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ABSTRACT 

We present the emission-line fluxes and kinematics of 48 representative elliptical and lenticu- 
lar galaxies obtained with our custom-built integral-field spectrograph SAURON operating on 
the Wilham Herschel Telescope. H/?, [O lll]AA4959,5007, and [N l]AA5 198,5200 emission 
lines were measured using a new procedure that simultaneously fits both the stellar spectrum 
and the emission lines. Using this technique we can detect emission lines down to an equiv- 
alent width of 0.1 A set by the current limitations in describing galaxy spectra with synthetic 
and real stellar templates, rather than by the quality of our spectra. Gas velocities and velocity 
dispersions are typically accurate to within 14 kms^^ and 20 kms^^, respectively, and at 
worse to within 25 kms^^ and 40 kms^^. The errors on the flux of the [O ill] and H/3 lines 
are on average 10% and 20%, respectively, and never exceed 30%. Emission is clearly de- 
tected in 75% of our sample galaxies, and comes in a variety of resolved spatial distributions 
and kinematic behaviours. A mild dependence on the Hubble type and galactic environment is 
observed, with higher detection rates in lenticular galaxies and field objects. More significant 
is that only 55% of the galaxies in the Virgo cluster exhibit clearly detected emission. The 
ionised-gas kinematics is rarely consistent with simple coplanar circular motions. However, 
the gas almost never displays completely irregular kinematics, generally showing coherent 
motions with smooth variations in angular momentum. In the majority of the cases the gas 
kinematics is decoupled from the stellar kinematics, and in half of the objects this decou- 
pling implies a recent acquisition of gaseous material. Over the entire sample however, the 
distribution of the mean misalignment values between stellar and gaseous angular momenta is 
inconsistent with a purely external origin. The distribution of kinematic misalignment values 
is found to be strongly dependent on the apparent flattening and the level of rotational support 
of galaxies, with flatter, fast rotating objects hosting preferentially co-rotating gaseous and 
stellar systems. In a third of the cases the distribution and kinematics of the gas underscores 
the presence of non-axisymmetric perturbations of the gravitational potential. Consistent with 
previous studies, the presence of dust features is always accompanied by gas emission while 
the converse is not always true. A considerable range of values for the [O lll]/H/3 ratio is 
found both across the sample and within single galaxies. Despite the limitations of this ratio 
as an emission-line diagnostic, this finding suggests either that a variety of mechanisms is 
responsible for the gas excitation in E and SO galaxies or that the metallicity of the interstellar 
material is quite heterogeneous. 
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1 INTRODUCTION 

Early-type galaxies were once considered uniform stellar systems 
with little gas, dust, and nuclear activity. A number of imaging and 
spectroscopic studies both from the ground and space have changed 
this view (see Goudfrooii 1999, for a review). Based on these sur- 
veys, we now know that early-type galaxies commonly contain dust 
in either organised or complex structures, whic h is almost always 
assoc i ated with optical nebular emission (e.g. , [Sadler & Gerhard 
1198.4 Ivan Dokkum & Franx 1995': IXran et alJ EoOll) . E arly-type 
galaxie s also show nuclear emission in 60% of the cases JHo et alJ 
Il997d) . 

Still, a number of issues remain open. What is the origin of 
the interstellar material in E/SO galaxies? Is it material lost by 
stars during their evolution o r does it have an external origin? It 
has long been demonstrated iFaber & Gallaghej[l973) that dur- 
ing their life stars would reinject more than enough material in 
the interstellar medium to explain the observed gas emission. Yet, 
the finding that the angular momentum of the gas or the orien- 
tation ofthedusti^veryoft^^ of the stars 
(e.g.. lBertola et alll992tlv'an Dokkum & Franxll995^ suggests an 
external origin. And what is its fate? Does it cool down to form 
stars or does it become hot. X-ray emitting gas? With the advent 
of new mm-wave detectors and the Chandra space telescope, both 
molecular gas and X-ray d etections have become more common 
in early-type galaxies (e.g.. I Younsll2002t lFabbiancll2003h . And fi- 
nally, what powers the observed nebular emission? Is it a central 
AGN? Is the warm (~ lO^K) gas ionised by the hot (~ lO^K) 
gas, through thermal conduction (the "evap oration flow" scenario, 
e.g- ISoarks et aJl98^lde Jong et alll990l) ? Is the gas ionised by 
stars, either young ('Shields" 1992^ or old (e.g., post-AGB stars, 
[di Serego Alishieri et al. 1990; Binette et al. 1994)? Or is the gas 
excited by shocks, as also pr oposed for low-ionisation nuclea r 
emission-line regions (LINERs. lPopita & Sutherlancil995lfl993) . 
If the extended emission observed in many early-type galaxi es is 
unlikely to be powered by nuclear activity iGoudfrooiill 1 999^ ■ all 
other ionising mechanisms are plausible. 

So far, the kinematics and ionisation of the gas in early- 
type galaxies have been studied mostly through long-slit obser- 
vations (e.g., Zeilinger et al. 1996; Caon et al. 2000; Phillios et al. 
1 198^ '. Ho et al. 1997b), while imaging surveys have investigated 
the distribution of the ionised gas and dust (e.g. , |Busoi^^1^^93|; 
iGoudfrooii 613011994 iMacchetto er^ll996l;Trran et alJl200ll) . 
Integral-field spectroscopic (IPS) data can combine both spatial and 
spectroscopic information, mapping the flux and kinematics of the 
ionised gas across large sections of nearby galaxies. Using two- 
dimensional measurements it is possible to accurately measure and 
compare the projected angular momentum of gas and stars, assess 
the regularity of the gas-velocity fields, and investigate the possible 
sources of ionisation for the gas. 

In this paper we present maps for the ionised-gas kinemat- 
ics and distribution within the effective radius of 48 representa- 
tive E and SO galaxies in both "cluster" and "field" environments, 
which were obtained with the SAURON integral-field spectrograph 
iBacon et aL.200L Paper I). These galaxies were observed in the 
course of the SAURON survey, a study of the structure of 72 repre- 
sentative nearby early-type galaxies and bulges. The objectives of 
the survey along with the definition an d properties of the sample 
are described in lde Zeeuw et ali ( 2002. Paper II). The stellar kine- 
matics for the 48 E an d SO galaxies in the survey are presented in 
lEmsellem et"an <2004 Paper III). 

This paper is organised as follows. In §2 we describe the 



extraction of the ionised-gas kinematics and fluxes, compare the 
SAURON measurements with published data, and set detection 
thresholds for the gas emission and the sensitivity of our survey. 
In §3 we present the maps for the ionised-gas distribution and kine- 
matics, discuss the incidence of emission, and describe the main 
features of the gas distribution and kinematics. In the same sec- 
tion we also describe the maps for the [O IIl]/H/9 line ratio and 
review systematic differences between the kinematics of the [O III] 
and H/3 lines. ij4|is devoted to the relation between gas and dust, 
while in ^we compare the kinematics of gas and stars in order to 
discuss the origin of the gas. We further discuss the observed gas 
phenomenology and draw our conclusions in §6. 



2 MEASURING THE IONISED-GAS KINEMATICS AND 
FLUX DISTRIBUTION 

Within the limited wavelength range of the SAURON observa- 
tions (4830-5330 A) there are three well-known emission lines 
we can expect to detect, the H/3A4861 Balmer line and the 
[OIII]AA4959,5007 and [N I] A A5 198,5200 forbidden-line dou- 
blets. The starting point for our emission-line measurements is de- 
scribed in Paper II. It consists of measuring the emission lines on 
residual spectra obtained by subtracting from the SAURON spectra 
a detailed description for the stellar spectrum, itself constructed for 
measuring the stellar kinematics. As described in Paper III and in 
Cappellari & Emsellem (2004), the latter process requires an op- 
timal combination of templates representative of the galaxy stellar 
population, excluding spectral regions that could be contaminated 
by emission lines. 

Further testing revealed that this procedure does not work 
well for all emission lines if there is insufficient information in the 
emission-free part of the spectrum to adequately constrain the stel- 
lar population content. Indeed when the wavelength range is limited 
it is possible that some of the emission lines we wish to measure lie 
very close to, or are coincident with, the most age- or metallicity- 
sensitive absorption features. By masking the regions potentially 
affected by emission, these important absorption features will be 
partially or even entirely excluded from the template-fitting pro- 
cess. This can lead to substantial biases in the resulting combination 
of templates, evident in the residual spectrum as spurious features 
in the masked regions that in turn contaminate the measurement 
of the emission lines. This is exactly the case for the SAURON ob- 
servations and particularly for the H/3 and [N I] doublet emission 
lines. The [O III] lines are less affected by this problem. 

2.1 The Method 

In this paper we extend the idea of Paper II, draw from the soft- 
ware of Cappellari & Emsellem (2004), and exploit the results of 
Paper III to develop a more refined procedure to measure the gas 
kinematics and fluxes without any spectral masking. The key in- 
gredient is to treat the emission lines as additional Gaussian tem- 
plates and, while iteratively searching for their best velocities and 
velocity dispersions, to solve linearly at each step for their ampli- 
tudes and the optimal combination of the stellar templates ^ , which 

^ The composition of the template library was impro ved w ith respect to 
the one used in Paper III. Three stars from the JoneJ il997h libraiy were 
exchanged with more suitable ones that allowed a better match to the spec- 
tra of large early-type galaxies, thus improving the emission-line measure- 
ments. The impact on the published stellar kinematics is negligible. 
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Figure 1. Examples of the difference between the emission-Hne measurements obtained when the lines are fitted to the residuals of the stellar-continuum fit, 
itself derived while masking regions contaminated by emission (top panels), or when both the stellar continuum and the emission lines are fitted simultaneously 
(lower panels). At the top of each panel the black histogram shows the galaxy spectrum, the red line the best fitting stellar template, and the green line the best 
fit when the emission lines are added. These are shown at the bottom of each panel by the blue line, along with a second histogram showing the difference 
between the galaxy specQum and the best stellar template. A constant has been added to both. The dashed vertical lines show the specfi'al regions that have 
been masked. Limiting the spectral range available to the template-fitting by masking introduces spurious features in the H/3 and [N I] regions that leads 
to overestimated amplitudes for these lines (top panels). The line widths are also overestimated, particularly when the [O III] lines are weak and no longer 
dominate the resulting gas kinematics (right). Both biases, in the Une amplitudes and kinematics, are overcome when the stellar continuum and emission fines 
are fitted together, and all the information in the spectra is used (lower panels). 



are convolved by the best stellar line-of-sight velocity distribution 
(LOSVD). In this way both the stellar continuum and the emission 
lines are fitted simultaneously. 

Here we adopt the same stellar kine matics and spa tial binning 
scheme of Paper III (from Icappellari & Copinll2003h and hence 
solve for the gas kinematics in exactly the same spectra. Further- 
more for each galaxy we adjust the continuum shape of the stellar 
templates using a multiplicative Legendre polynomial of the same 
order (typically 6) as needed in Paper III. We adopt a multiplicative 
polynomial adjustment to ensure that no extra dilution of the ab- 
sorption features is introduced. As mentioned before, we assume a 
Gaussian LOSVD for the gas clouds. In the case of doublets, each 
component has the same mean velocity and width, and their relative 
strength is fixed by the ratio of the corresponding transition prob- 
abilities. Figure Q shows for two specific cases the advantages of 
this new method, when all emission-lines are fitted with the same 
kinematics. Across our sample, template-mismatch features in the 
masked regions lead to overestimated H/3 fluxes in average by 10% 
and by up to 40% and, if the [O III] lines are weak, also to system- 
atically overestimated line-widths by ~ 10%. 

In principle, we could have used the new method imposing 
the same kinematics on lines emitted from different atomic species 
or while fitting each of them independently. In practise, however, 
neither the H/3 nor the [N I] lines could be always measured confi- 
dently without first constraining their kinematics, as contamination 
due to template mismatch can still be important even using our new 
approach. 



The main difficulty with the H/3 measurement is the presence 
of a number of metal features, mainly from chromium and iron, 
around 4870A. If we allow the position and width of the H/3 emis- 
sion to vary, when the stellar templates cannot match the strength of 
this absorption feature, the overall spectral shape between 4855A 
and 4875A can be described using metal-poor templates and a spu- 
rious H/3 line placed roughly half-way between the H/3 and Cr-Fe 
absorption features. Figure |2| shows an example of this problem. 
The strength of this emission artifact can be significant and produce 
a spurious detection. Fortunately, because the spurious H/3 emis- 
sion always falls at the same wavelength region in the rest frame, 
the mean velocity of this line appears shifted with respect to the 
stellar velocity by an approximately constant amount. Hence it is 
possible to recognise this problem across the field by comparing 
H/3 velocity maps to the stellar velocity maps (Figure|2j- 

The lines of the [N I] doublet are normally quite weak so that 
it is almost never possible to constrain their kinematics indepen- 
dently. To complicate matters further, the [N I] doublet sits close 
to the continuum region that is generally the worst matched by 
the templates, because of an enhanceme nt in the magnesi um over 
iron ratio that is neither included in the IVazdekij i 19991) models 
(which form the bulk of our template library) nor observed in stars 
in our Solar neighbourhood. The impact of template-mismatch on 
the measurement of the [N I] lines is therefore more difficult to es- 
timate than in the case of the H/3 emission. 

To measure the fluxes of the H/3 and [N I] lines we therefore 
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Figure 2. Example of template-mismatch contamination in the H/3 measurement, in the SO galaxy NGC 3489. (a) Spectral region around H/3. The soHd black 
histograms show the galaxy spectrum, the green lines the best fitting stellar template, and the red lines the best fit when the H/3 emission, indicated below by 
the solid lines, is added. When the velocity of H/3 emission is free to vary (left) the two main absorption features are best fit by a combination of templates and 
a spurious H/3 emission feature that contributes to the 'bump' in the spectrum between the two absorptions. When the velocity of H/3 emission is constrained 
to that of [O III] (light), the best fit includes a weaker, naiTower emission feature closer to the H/3 absorption, (b) Velocity maps. As the unconstrained H/3 
emission is a spurious feature roughly always at the same relative position with respect to the H/3 absorption feature, approximately 2/3 of a pixel (60km/s) to 
the red, the corresponding velocity field (left) resembles the stellar velocity field when a constant is added to it (center, note the curved zero-velocity region in 
green). The right panel shows the [O III] emission velocities, which trace the real gas velocities. 



constrained their kinematics to those of the [O III] lines, which was 
obtained first. 

Using the spatial binning, spectral library, and stellar kinemat- 
ics of Paper III, we obtain the emission-line fluxes and kinemat- 
ics in each (binned) spectrum of our sample galaxies by following 
these exact steps: 

a) IVIask all spectral regions within ±300kms~^ of the location 
of the H/3 and [Nl] lines at the galactic systemic velocity Vsys, 
taken from Paper III. This covers the typical range of emission-line 
velocities and widths found in a preliminary analysis of our sample. 

b) Convolve all stellar templates in our library with the corre- 
sponding best stellar LOS VD from Paper III. 

c) Solve for the best amplitude y4gas, mean velocity 1/gas and in- 
trinsic velocity dispersion dgas of the [O III] lines, while also opti- 
mising the continuum shape of the templates using a multiplicative 
polynomial adjustment. Vsys is used as initial velocity guess for 

d) Remove the mask and find the best Agas for all the lines, while 
imposing on them the [O III] kinematics. 

IVIore specifically, at step c) the best Vgas and cTgas and the 
best coefficients of the multiplicative polynomial correction are 
found through a Levenberg-Marquardt least-squares minimisation. 
At each iteration we construct a Gaussian template for each emis- 
sion line with the current position and width (accounting for the 
SAURON spectral resolution), and with unit amplitude. For the 
[O III] doublet, the template is formed b y two Gaussians, with am - 
pUtude of 0.33 for the [O III]A4959 line jStorev & ZeiDpeij2000l) . 
After multiplying the convolved stellar templates by the current 
polynomial correction, we fit for the best linear combination of both 
stellar and emission-line templates (with positive weights), exclud- 
ing the regions potentially affected by H/3 and [N I] emission us- 
ing the mask built in a). The weights assigned to the emission-line 
templates provide the best emission-line amplitudes Agas at each 
iteration, and eventually, the final best values. 

Step d) is similar to c), although only the coefficients of the 
polynomial adjustment are solved for non-linearly, while the con- 
tribution of the stellar and emission-line templates is still optimised 
at each iteration. Gaussian templates are constructed also for the 
H/3 and [N I] lines, and the entire spectrum is used in the fit. 

In order to follow the most general method described at the 
beginning of this section, only steps b) and c) are needed, without 
masking any spectral region and by searching for the gas kinemat- 
ics of all lines at the same time. 



2.2 Constructing the Emission-Line IMaps 

In our sample galaxies the ionised-gas emission is neither uni- 
formly distributed nor always strong enough to be detectable. It is 
therefore crucial to understand the level to which we are confident 
of detecting gas emission. In Appendix A we present a number of 
experiments specifically designed to address this issue, and here we 
summarise the results obtained there. 

The accuracy with which the position and width of an 
emission-line can be recovered depends on how much the line pro- 
trudes above the noise in the stellar spectrum. We call this quantity 
the line amplitude-to-noise ratio, hereafter A/N. The accuracy in 
recovering the amplitude of a line, on the other hand, scales only 
with the noise level in the spectrum. The ability to estimate the am- 
plitude of the lines is also the dominant factor in the error budget of 
the line fluxes. Hence, in the limit of purely statistical fluctuations 
the line fluxes are subject to larger errors in spectra of higher qual- 
ity, although the equivalent width of the lines is better estimated at 
these regimes. Better spectra also allow a more accurate description 
of the stellar continuum, which is equally crucial to the emission- 
line measurements. 

Indeed, since emission lines are measured while simultane- 
ously fitting the stellar spectrum using a template library, any sys- 
tematic mismatch between the templates and the galaxy stellar pop- 
ulation will constitute a further source of error. It is therefore im- 
portant to include such deviations when estimating the noise level 
against which the emission-l ine ampl itudes are compared. We use a 
robust biweight estimator ( Hoaglin et a l. 1983) to measure the scat- 
ter in the residuals of the fit to the integrated stellar spectrum, as an 
estimate of both statistical fluctuations and systematic deviations. 
We will refer to this as to the "residual noise". 

The simulations of Appendix A also show that the formal un- 
certainties returned by our emission-line fitting procedure correctly 
estimate the accuracy with which the input parameters are recov- 
ered. At low A/N, however, our measurements become dominated 
by systematic, rather than random, effects and the formal uncertain- 
ties cannot account for the observed biases. The minimum A/N 
values below which these problems arise naturally sets our detec- 
tion threshold. 

According to our emission-line measurement scheme, we 
must first assess the presence of the [O III] emission, without which 
no other emission line can be detected. The simulations show that 
our ability to estimate the velocity dispersion of the [O III] lines 
steadily deteriorates for A/N ^ 2.5, while the line fluxes start 



The SA UR ON project - V 5 



to be overestimated by more than 10% for A/N ^ 5. As these 
experiments cannot fully account for the limitations of our tem- 
plate library, we conservatively settled on an overall A/N threshold 
of 4. This limit was chosen also considering that in galaxies with 
larger velocity dispersions the impact of template-mismatch can be 
more important and cause the line widths to be overestimated, as 
shown by our experiment based on independent fits to the two lines 
in the [Olll] doublet f ^A5r For an intrinsic cfgas = 50km s"^ 



and A/N — 4 the simulations show that the errors on l^as and 
(jgas are typically 25 kms^^ and 40 kms~^, respectively, while 
errors on fluxes are 30%. The large error on the intrinsic (jgas is 
not due to a poor fit but arises because the instrumental resolution 
is much larger than the adopted value of 50 kms~^. In fact, the 
observed width of the lines is matched to within 17 km s~^. Mea- 
surements for intrinsically broader lines will be more accurate for 
a given A/N. 

When [O III] lines are detected, H/3 emission is also likely to 
be present. The possibility that the H/3 emission could be dom- 
inated by gas with different kinematics than that contributing to 
the [O III] emission will not change our ability to estimate the H/3 
fluxes. Indeed, the accuracy in recovering fluxes is not dramati- 
cally affected by how well the position and width of the lines are 
measured. This is not only observed in our simulations, but also 
in the few galaxies where the [O III] and H/3 lines could be fitted 
independently (' il3.5> . When we impose the [O III] kinematics on 
the H/3 line, our simulations show that, for A/N ^ 2.5, the esti- 
mates of the H/3 flux are unbiased. At this level, they are accurate 
to within ~30%. Unbiased H/3 fluxes can be measured with lower 
A/N because the H/3 fit involves a lower number of parameters 
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Figure 3. SAURON vs. Palomar. [O III]/H/3 line ratio ( left) and intrinsic 
width (right) of the [N II] lines from the nuclear spectra of lHo et al]ll997hll 
compared to [O IIIJ/H/3 ratio and intrinsic width of the [O III] lines mea- 
sured with SAURON in the same apertures. The lower panels show the two 
sets of data directly against each other, while the upper panels show their 
ratio. Larger symbols correspond to less accurate measurements from Ho et 
al., while the aiTows con'espond to upper and lower limits from the same 
work. The dashed lines show the identity relation, while the dotted lines in 
the top panels define the region where the SAURON and Palomar measure- 
ments agree to within 20%. 



fects the measured H/3 fluxes even in the framework of ideal sim- 
ulation, we conservatively settled on a A/N threshold of 3 for this 
line. 

Finally, we come to the harder problem of assessing the pres- 
ence of the [N I] doublet. These lines are the most affected by our 
limited ability to describe the stellar populations of our sample 
galaxies, a problem that we cannot properly simulate. The most 
likely source of template mismatch in this spectral region is the 
absence of templates with super-solar abundance ratios. A more 
robust detection threshold can be established by comparing the 
[Nl]A5200 amplitude to the residual-noise level measured more 
specifically in the Mg region, A'^Mg. Furthermore, we do not ex- 
pect to observe [N I] lines without strong H/3 and [O III] emission. 
Hence, we consider the detection of the [N I] doublet to be reliable 
only if both the H/3 and [O III] lines have already been detected, 
and if the [N I] lines satisfy the quite conservative detection re- 
quirement of A/A^Mg ^ 4. 

For our survey we will therefore show maps for the flux, 
equivalent widths, velocity, and velocity dispersion of the [O III] 
emission-lines for which A/N 4. Maps for the flux and equiv- 
alent width of the H/3 lines and for the [O IIl]/H/3 ratio will show 
regions where additionally A/N ^ 3 for the H/3 line. We will not 
show maps for the [N I] lines, as these are detected only in the cen- 
tral regions of 13 objects. Instead, the galaxies with [N I] emission 
are listed in TableQ 

To conclude, we stress that the detection thresholds adopted 
here will never perfectly exclude all spurious measurements or 
guarantee that all regions with real emission appear the maps. For 
completeness, the public data release will contain all emission-line 
measurements with associated errors, for all spectra corresponding 
to single lenses or larger spatial bins shown in the maps, with a flag 
marking the measurements we deem unreliable. 



2.3 Sensitivity Limits 

The detection thresholds adopted in the previous section set the 
sensitivity of our survey. The equivalent width {EW) of the weak- 
est [O III] and H/3 lines that we detect are around O.I A and O.OTA 
for the [O III] A5007 and H/3 lines, respectively. These limits can be 
understood as follows. If S is the continuum level in the spectra, F 
the flux in the lines, and a their typical observed widths, then an 
A /N threshold can be translated into a limiting EW considering 
that for a Gaussian line: 



EW ■ 



F 



A X 



S 



2na A/N X 

or EW = ' 



S/N 



(1) 



Hence, for emission lines with a typical intrinsic broadening of 
cTgas = 50 kms^^, which due to instrumental broadening appears 
as a = 120 kms"\ or 2A at 5007A, the EW of a barely detected 
[O III]A5007 line (A/N = 4) would be 0.2A for S/N = 100. 

In the nuclear region of our sample galaxies the SAURON spec- 
tra are of extremely high quality, with nominal 5'/A*' up to 500 per 
pixel, so that very weak lines should be detected. However, what 
matters for the detection of emission lines is the "residual noise", 
not just the statistical fluctuations in the stellar spectrum. When 
the "residual noise" is compared to the continuum level, the corre- 
sponding 5'/A'^ ratios reach values only up to ~ 200, dominated by 
template mismatch. With this upper limit for S/N the EW limits 
become 0.1 A for [O III]A5007 and 0.07A for H/3, as found. At a 
lower surface brightness level, however, the statistical fluctuations 
still dominate the "residual noise", and the signal-to-noise ratio is 
close to nominal. In Paper III we adopted a binning scheme with 
5'/A'^ ^ 60. Therefore our sensitivity is never worse than 0.3A and 
O.2A for [O III]A5007 and H/3, respectively. 
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2.4 Comparison with Published Emission-Line 
Measurements 

In Paper II we showed for NGC 5813 that the SAURON ionised- 
gas kinematics is consistent with published data. An ideal source 
for a more general comparison is the Palomar spectroscopic sur- 
vey of Ho, Filippenko, & Sargent (1995, 1997a), who observed 
37 of our 48 sample galaxies and also carefully subtracted the 
stellar light prior to the emission-line measurements. We com- 
pare the [O IIl]/H/3 emission-line ratios as well as the width of 
the forbidden emission, in this case r elating the SAUR ON [O III] 
lines width to the [N II] A6583 width of lHo et alJ <1997ah . For con- 
sistency with the Palomar data, we analysed nuclear spectra ex- 
tracted from our SAURON cubes within central apertures that match 
the size (2" x 4") and orientation of the Palomar long-slit ob- 
servations. For these spectra we first derived the stellar kinemat- 
ics as in Paper III, using the penalised pixel fitting algorithm of 
ICappellari & EmsellemI i2004h . and then followed our procedure 
to measure the gas emission. 

The result of this comparison is shown in Figure |3| for the 
21 objects with measured emission-lines in the Palomar survey. 
Considering the number of possible systematic factors that could 
enter this exercise (e.g., different observing conditions, data qual- 
ity and reduction, starlight subtraction), the agreement between the 
SAURON and Palomar measurements is satisfactory, as these are on 
average consistent to within 20% (see Figure|3). In particular, for 
the [O III] over H/3 ratios, the quality of the Palomar blue spectra 
is significantly worse than that of the SAURON data, with nominal 
S'/A'^ ^ 20 — 30 per pixel as opposed to at least 500 for our nuclear 
extractions. On the other hand, the higher quality of the Palomar 
red spectra {S/N ~ 30 — 80) is consistent with the better match 
between the widths of the [O ill] and [N II] lines. 

Forbidden lines can have different widths depending on 
their respective criti cal densities ( Filiooenko & Haloern 1984; 
Ide Robertis & Osterbr ock 1986). The fact that in the majority of 
the SAURON galaxies the width of the [N II] lines is compara- 
ble to that of the [O ill] lines could be explained if emission pre- 
dominantly originates in low-density reservoirs, as suggested by 
[Ho et al . ( 1997b) in the case of [Nil] and [S II], which also often 
display similar widths. On the other hand, the origin of the four 
outliers in Figure |3| with significantly broader [Olll] lines can be 
explained considering that, according to the positive correlation be- 
tween line width and critical density, [O III] should be broader than 
[Nil]. 

Of the 16 objects that Ho et al. identified as emission-line free, 
9 show clear emission in the SAURON central apertures. Apparently 
low-luminosity nuclear activity in early-type galaxies is even more 
common than already established. 



3 lONISED-GAS DISTRIBUTION AND KINEMATICS 

Fi gure I4al4bl shows maps for the flux and equivalent widths of the 
[O III]A5007 and H/3 emission lines, for the velocity and intrinsic 
velocity dispersion of the [O III] fines, and for the [O IIl]A5007/H/3 
ratio of all 48 E and SO galaxies in our sample. The maps were 
constructed according to the detection thresholds set in i|2.2l The 
SAURON observations were not always carried out in photometric 
conditions. Therefore the derived fluxes should be considered only 
as approximate values. 

Appendix lElcontains a description of the gas distribution and 
kinematics in each of our sample galaxies, along with comments on 



the dust distribution, the connection to the stellar kinematics, the 
[O IIl]/H/3 ratios, and references to previous narrow-band imaging 
and long-slit spectroscopic work. In the following we quantify the 
incidence of ionised-gas emission and describe the general prop- 
erties of the gas distribution, the morphology of the gas velocity 
and velocity dispersion field, and the range of [O IIl]/H/3 ratios ob- 
served in Figure l4al4bl for galaxies with clear emission-line detec- 
tions. We will also briefly investigate the presence of systematic 
differences between the kinematics of the [O ill] and H/3 lines. 

3.1 Incidence of lonised-Gas Emission 

The maps reveal that ionised-gas emission is very common in 
early-type galaxies, and that it comes with a variety of spatial dis- 
tributions, degrees of regularity in the observed kinematics, and 
large variations in the [O IIl]A5007/H/3 ratio. [O III] and H/3 emis- 
sion is clearly detected in 36/48 galaxies in our sample (75%), 
while in a further 7 objects we found either weak central [O III] 
lines only (NGC 4382, NGC 5845) or patchy traces of emission 
(NGC 4270, NGC 4458, NGC 4473, NGC 4621, NGC 4660). 
Whether strong, weak or patchy, such emission is always spa- 
tially resolved. Just 5 galaxies do not show any significant emission 
(NGC 821, NGC 2695, NGC 4387, NGC 4564, NGC 5308). [N I] 
lines are found in 13 galaxies. TableQlists the basic properties of 
the E/SO SAURON sample and identifies galaxies in which H/3 and 
[O III] emission was detected, and where [N I] lines were found. 

Among the clear detections, the incidence of emission lines 
is higher in lenticular galaxies, where emission is found in 20/24 
objects (83%), compared to 16/24 (66%) for ellipticals. The depen- 
dence on the galactic environment appears similarly marginal, with 
20/24 field galaxies showing emission compared to 16/24 in clus- 
ters, where the definition of "field" and "cluster" is as in Paper II. 
However, the fraction of galaxies with clearly detected emission in 
the Virgo cluster drops to only 55% (10/18), with just 3/9 ellipticals 
showing the presence of ionised gas. These 3 galaxies are also the 
brightest ellipticals that we observed in Virgo. More generally the 
incidence of emission does not seem to depend on the galaxy lu- 
minosity. The incidence of gas emission is the same for barred and 
unbarred SO galaxies. 

Table Q also lists the total flux of the H/3 emission in each 
galaxy with clearly detected emission. This includes the contribu- 
tion of regions where only the [O III] lines were detected, assuming 
a constant line ratio of [O IIl]A5007/H/3 = 3. Assuming Case B 
recombination, a temperature of T = 10^ K, and an electron den- 
sity of rie — 10^ cm^'', we als o rep ort total Ha luminosities and 
ionised-gas masses (following lKini ll989). In computing the Ha 
luminosity we adopted the distance moduli of Paper II, and used 
the th eoretical va lue of 2.86 for the flux ratio of the Ha and H/3 
lines (Osterbrock 1989), thus ignoring the impact of dust absorp- 
tion which would lead to higher Ha luminosities and larger gas 
masses. We caution against using such flux and mass measurements 
for quantitative applications. 

3.2 lonised-Gas Distribution 

The distribution of the ionised-gas is traced by the flux maps in 
Figure|4a| Figure |4b| also includes maps for the EW of the lines, 
which highlight structures in the emission-line distribution that are 
otherwise hidden in the flux maps (e.g., NGC 2974, NGC 3414). 
This is because the flux distribution of the lines mostly follows that 
of the stellar continuum, with a wide dynamic range. Small fluctu- 
ations are therefore harder to recognise in the flux maps, whereas 
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the EW -maps have the variation of starlight divided out. Although 
EW structures are not directly related to the flux distribution of 
the lines, they show specific emission regions such as rings or spi- 
ral arms. The EW -maps also provide a picture of the robustness 
of our emission-line measurements and illustrate how close they 
come to the detection limit. The EW is indeed closely related to 
the A/N ratio (Eq.0, which sets our detection thresholds and de- 
fines the accuracy of the kinematic measurements. The EW-maps 
should be used only as supplementary tools to the flux maps when 
investigating the gas distribution. Overall, by complementing the 
flux maps with the EW-maps, a better picture of the ionised-gas 
distribution emerges. 

The simplest case is that of an extended distribution consistent 
with a disk of gas. NGC 524, NGC 4459, NGC 4526, NGC 5838 
belong to this category, as do the polar-ring galaxies NGC 2685 and 
NGC 2768 where the gas distribution is perpendicular to the stel- 
lar body. Additional spiral or ring features are found in NGC 2974, 
NGC 3414, NGC 4550, and possibly in NGC 3608. In some other 
galaxies (NGC 5198, NGC 5846, NGC 5982) an elongated struc- 
ture extends from a central disk. A very distinct class of objects is 
formed by galaxies where a loosely-wound spiral feature, similar 
to an integral sign, is superimposed on a more diffuse component 
(e.g., NGC 474, NGC 3377, NGC 4262, NGC 4278, NGC 4546). 
Among more complex distributions are NGC 3156 and NGC 3489, 
which exhibit a central ring and filamentary structures, and objects 
characterised by outer emission-line regions that are strongly mis- 
aligned with the main body of the galaxy (NGC 1023, NGC 2549, 
NGC 7332, NGC 7457). Finally, there are giant ellipticals where 
the gas distribution is mainly confined to lanes across the galaxy 
(NGC 4374, NGC 5813) or in filaments (NGC 4486). 

We note that planetary nebulae (PNe) could be responsible 
for many of the isolated patches of [O III] emission observed in 
our sample galaxies. For instance, the typical observed fluxes of 
the PNe in NGC 3379 ( ~ 10"^'^ erg s~ for a narrow-band 

magnitude msoo? = 26, ICiardullo etalll989h would correspond 
to equivalent widths well above our sensitivity limit, except in the 
very central region (within 5") and for the largest outer bins. The 
detection of PNe in NGC 3379 is discussed in AppendixlBl 



3.3 lonised-Gas Kinematics 

3.3.1 Velocities 

A number of galaxies clearly show very regular gas veloc- 
ity fields consistent with circular gas motions in a disk (e.g., 
NGC 524, NGC 4459, NGC 4526). Yet, the majority of the ob- 
jects appear to deviate from this simple situation, showing grad- 
ual twisting of the overall velocity field (NGC 2974, NGC 3377, 
NGC 3414, NGC 4278), more complex leading or trailing features 
(e.g., NGC 2768, NGC 4550), and sometimes both (NGC 3489, 
NGC 4374). More mundanely, heavy spatial binning or weak emis- 
sion leads to more noisy velocity fields (e.g., NGC 4150). In some 
objects the radial variation of the gas angular momentum is more 
abrupt. NGC 474 is the clearest example of this behaviour, but also 
NGC 2549, NGC 7332, NGC 7457, and possibly NGC 1023, show 
kinematics that are distinct in the inner and outer regions. Overall 
across our sample the gas kinematics are never irregular, the only 
exception being NGC 4486. 

We note that regular motions occurs preferentially in objects 
with a regular disk distribution, while a more complex kinematic 
behaviour follows the presence of additional features in the gas 
distribution, such as integral-sign and spiral structures. In partic- 



ular, all objects with relatively strong, misaligned emission in their 
outer parts (e.g., NGC 2549), also show strongly decoupled exter- 
nal and internal gas motions, with the innermost system being more 
likely settled onto the equatorial plane. 

3.3.2 Velocity Dispersions 

The accuracy of the (jgas measurements is more sensitive to the 
impact of weak emission than for the V^as measurements. This is 
exhibited as an increasing fluctuation in o-gas in regions with barely 
detected [O III] emission (e.g., NGC 4374). At low A/N the width 
of the lines also tends to be overestimated, due to template mis- 
match and correlation in the random fluctuations in the spectra (see 
SjASland JAS). 

The (Tgas-maps appear more uniform than the EW and Igas 
maps, and are characterised by little variation of dgas over large 
regions and by the presence of a central gradient. Both the typi- 
cal width of the lines in the outer regions and the magnitude of 
their central increase can vary considerably. This is evident when 
the (Tgas -values are plotted against their distance from the centre 
(Figure |4bJ. In some cases the intrinsic dispersion of the lines re- 
mains as high as 100 kms^^ in the outer parts (e.g., NGC 4278, 
NGC 4546) while in other falls to 50 kms"^ (e.g., NGC 3489, 
NGC 3156). Central values of (jgas reach ~200 kms~^ (e.g., 
NGC 2974, NGC 3414). In other cases f7gas remains constant 
(NGC 3489). The morphology of the central peaks can vary, from 
shai-p (e.g., NGC 3377) to very extended (e.g., NGC 2974), elon- 
gated (e.g., NGC 2768, NGC 3414), and even asymmetric (e.g., 
NGC 4262, NGC 4278, NGC 4546). Finally, some objects show 
pecuhai- (jgas -features (e.g., NGC 474 and NGC 2685). 

Except for asymmetric (jgas-peaks that occur in objects with 
integral-sign _EW^-patterns, the most significant characteristics of 
the (Jgas -maps do not appear to correlate with the morphology of 
the gas distribution and velocity field. In the central regions, how- 
ever, the observed velocities are important to understand whether 
the (Jgas -gradients are due to unresolved rotation, an AGN, or a 
genuine increase of (Jgas. Detailed models are needed to address 
this issue, but we expect the effects of purely unresolved rotation 
to be rather limited and confined to the seeing-dominated regions. 
The strongest and most extended (Jgas -gra dients therefore sugges t 
an intrinsic rise of (jgas (e.g., NGC 2974. iKrainovic et alJl200 j) . 
On the other hand, the case for unresolved rotation is particularly 
strong for sharp peaks that are elongated along the direction of the 
zero-velocity curve (e.g., NGC 2768, NGC 3414). 

We note that (jgas always exceeds the expected value from 
thermal broadening (~ 10km ^ ) and is generally smaller than the 
stellar velocity dispersion cr,. In some cases, however, (jgas ~ (J*, 
either only in the central regions (e.g., NGC 2974, NGC 3414, 
NGC 5813) or over most of the field (e.g., NGC 3156). This sug- 
gests the presence of additional turbulence in the gas clouds or that 
the latter do not follow perfectly c ircular orbits. A stro nger dynam- 
ical support for the gas motions iBertola et ai]ll995h may be re- 
quired to explain cases where (jgas ~ cr* . 
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Figure 4a. Maps of the ionised-gas distribution and kinematics of tiie 48 E and SO galaxies in the SAURON representative sample. The SAURON spectra 
have been spatially binned to a minimum S/N of 60 in the stellar continuum, consistent with Paper III. All maps are plotted on the same angular scale, in 
arcseconds. The aiTow and its associated dash at the top of each column mark the North and East directions, respectively. From top to bottom: i) and ii) flux 
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Figure 4b. From top to bottom: i) and ii) equivalent width of the [O III]A5007 and H/3 emission lines, in A and in a logarithmic scale, iii) stellar mean velocity 
in km from Paper III, iv) unshaip-masked images obtained from HST observations or SAURON reconstructed intensity maps, and v) radial profiles for the 
misalignment between the kinematics of gas and stars (up), for the velocity dispersion of the [O III] lines (middle) and for the [O III]A5007/H/3 ratio (down). 
The SAURON maps are as in Figure l4al The grey boxes in the top two maps indicate the field of the HST images. 
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3.4 [O iii]/H/3 ratios 

The maps for the relative strength of the [O III]A5007 and H/3 hnes 
in Figure|4a|can be used to identify regions where emission could 
be powered by young stars and to trace variations of the ionisation 
mechanism within single galaxies. Low [O IIl]/H/3 ratios tend to 
characterise star forming H Il-regions, whereas other mechanisms 
are in general responsible for emission with high [O III]/H/? ra- 
tios. Only ~20% of the emission-line nuclei with [O IIl]/H/3> 1 
are classified as Hll nuclei, while this class of objects represent 
~85% of the nuclei with [O IIl]/H/3s; 1 (Hoetal. 1997c). On 
the other hand, [O III] lines stronger than H/3 can arise also in 
Hll-regions if the metallicity of the gas is sufficiently low (e.g., 
[Veilleux & OsteAroc3[l98J. The metallicity of the interstellar 
medium can vary between different objects (e.g., if the gas has an 
external origin) but it is unlikely to change abruptly across different 
regions of a galaxy. Large fluctuations of the [O IIl]/H/3 ratio within 
a galaxy suggest a variation of the ionising mechanism rather than 
a change in the gas metallicity alone. 

The great variety of [O IIl]/H/3 values found both across our 
sample and within single objects (Figurel4al therefore suggests ei- 
ther that different ionising sources could be at work in early-type 
galaxies or that the metallicity of the interstellar medium is very 
heterogeneous from galaxy to galaxy. The radial profiles of Fi d4bl 
show that more than half of the objects have average [O IIl]/H/3 val- 
ues between 1 and 3. Many galaxies (e.g., NGC 3377, NGC 3489) 
display emission with [O IIl]/H/3> 3 , but such a strong level of 
excitation (typical of Seyfert nuclei) never dominates the entire 
field (except for NGC 7332). On the other hand, [O IIl]/H/3 ra- 
tios ^ 1 characterise most of the emission observed in NGC 524, 
NGC 3032, NGC 4459 and NGC 4526. The radial profiles for the 
[O IIl]/H/3 ratios also emphasise the presence of central gradients. 
Towards the centre, with few exceptions, [O IIl]/H/3 always ranges 
between 1 and 3. 

The ionisation structure we observe links to the gas distribu- 
tion and kinematics in two classes of objects in particular. The first 
group includes the four galaxies with the lowest [O IIl]/H/3 ratio 
in our sample: NGC 524, NGC 3032, NGC 4459 and NGC 4526. 
These objects have clear disk-like gas distribution and kinemat- 
ics, although in NGC 3032 severe binning makes it hard to judge 
the regularity of the velocity field. NGC 3032, NGC 4459, and 
NGC 4526 also share the same radial trend for the [O IIl]/H/3 ratio, 
with a minimum in a circumnuclear region and stronger [O IIl]/H/3 
ratios toward the centre and the edge of the gas disk. Emission in 
NGC 524 is too weak to confirm this behaviour. The second class is 
composed of NGC 4262, NGC 4278, and NGC 4546, which show 
[O IIl]/H/3-maps that are not symmetric around their centre despite 
the corresponding [O III] and H/9 distributions being fairly simi- 
lar. These objects also display an integral-sign feature in their gas 
distribution. By contrast, NGC 3377, another object with this sig- 
nature, has a fairly symmetric [O IIl]/H/3 distribution. 

We note that many of the brightest elliptical galaxies in our 
sample show rather uniform [O IIl]/H/3-maps for relatively inter- 
mediate values of [O IIl]/H/3~l-2. 



3.5 H/3 and [O ill] kinematics 

As clouds ionised by different sources need not share the same 
kinematics, here we seek evidence for systematic differences be- 
tween the kinematics of clouds that emit more efficiently H/3 pho- 
tons than [O III]. We measured the H/3 and [O III] kinematics inde- 
pendently in the galaxies with the strongest emission, and checked 
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Figure 5. Faster rotating H/3 kinematics in galaxies with perfect dusty disks. 
Two examples of galaxies with very regular gas kinematics and circularly 
symmetric dust lanes where the [O III] and H/3 kinematics could be derived 
independently. The velocity maps for the [O III] (top) and H/3 lines (mid- 
dle) are as in Fi d4al and show regions where both [O III] and H/3 lines were 
detected following our standard approach. Lighter and darker colors on the 
receding and approaching side, respectively, of the H/3 velocity maps, illus- 
trate how the H/3 kinematics show faster rotation velocities than the [O III] 
kinematics. The grey boxes in the velocity maps indicate the field of the 
HST images (bottom), showing the dust-lane morphology (see also Fi j4bl . 



if the H/3 measurements were subject to the systematics described 
in §2.1 by comparing the H/3 velocity field with the stellar one 
(see Figure |2j. The [N I] lines were still fitted using the [O III] 
kinematics. Independent H/3 and [O III] kinematics could be mea- 
sured over most of the region where emission is observed only 
in 10 galaxies. The most significant result from this experiment, 
shown in full in Appendix |C| is that in objects with very regular 
gas distribution and kinematics, such as NGC 4459 and NGC 4526, 
the H/3 velocities are higher than the [O III] velocities, by up to 
100 kin s^^(Figure|5}. This suggests that H/3 is a better tracer than 
[O III] for the circular velocity in the equatorial plane. 

In other cases the differences between the velocities are more 
limited (^ 30km s~^) and show less clear patterns, except for the 
central regions of NGC 3414 and NGC 4278 where a faster [O III] 
component is observed. In general the width of the H/3 lines tend 
to be smaller than that of the [O III] lines, in particular towards the 
centre. 

We note that even when the difference between the H/3 and 
[O III] kinematics is extreme (100 kms^^), the flux of the H/3 lines 
fitted independently and of those measured with the [O III] kine- 
matics differ by less than 20%. In terms of equivalent width the 
constrained fit can lead, in the most extreme cases, to an overesti- 
mation of H/3 by up to 0.2A. 
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Table 1. lonised-gas emission tlie 48 E/SO SAURON galaxies. 
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Notes: (1) NGC number (2) Hubble type (RC3: lde Vaucouleurs et alll99lll . (3)— (5) Galactic Environment, distance modulus Am in mag, and absolute 
_B-band magnitude, from Paper 11. (6) Ellipticity £25 of the 25 _B-band mag arcsec"^ isophote (LEDA). (7) Presence of H/3 or [O III] AA4959,5007 emission. 
(8) Presence of [N I]AA5198,5200 lines. (9) Presence of dust features in the HST images. (10) - (12) Estimated (see text) total H/3 flux. Ha luminosity and 
mass of the ionised gas in erg s~^cm~^, erg s~^, and M©, respectively, for objects with clearly detected emission. 



4 DUST AND GAS 

Dust in early- type galaxie s is alm ost always asso ciated with gas 
emission (e.g. jTomita et a l. 2000; T ran et all200ll) . With SAURON 
we can correlate the presence and spatial morphology of the dust 
not only with the ionised-gas distribution but also with its two- 
dimensional kinematics. 



To highlight the presence of dust we constructed unsharp- 
masked maps using both the SAURON reconstructed intensity maps 
and arciiival HST images, the latter being available for 45/48 ob- 
jects. This was done by dividing each image by a smoothed ver- 
sion of itself, using a Gaussian kernel with a FWHM of (/.'S and 
(X'S for the HST and SAURON images, respectively. Figure I4bl 
shows for each galaxy the HST or SAURON unsharp-masked im- 
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age. The unsharp-masked images reveal different dust morpholo- 
gies, such as perfect disks (NGC 524, NGC 3032, NGC 3379, 
NGC 4459, NGC 4526, NGC 5838), less defined coplanar distri- 
butions, (e.g., NGC 2974, NGC 3489, NGC 4550), well defined 
lanes (e.g., NGC 2685, NGC 3377, NGC 4374) or more complex 
filamentary (e.g, NGC 4486, NGC 5813, NGC 5846) and patchy 
(e.g., NGC 4552) structures. Whenever dust is present we also find 
emission, although very weak in NGC 5845. The converse does not 
always hold, which is expected as dust is harder to detect in close to 
face-on configurations. Still, it is puzzling to see that in some cases 
(e.g., NGC 5198, NGC 5982) dust is not detected despite the ob- 
served high gas velocities, which are hard to reconcile with motions 
in a face-on disk. 

Consistent with previous studies (e.g., iGoudfrooii et all 1994) . 
the dust generally follows the ionised-gas distribution, even if it is 
not always possible to connect features in the flux and iJVF-maps 
with ones in the unsharp-masked images. In some objects, however, 
specific H/3-emitting regions have dusty counterparts. Perfect ex- 
amples are the circumnuclear H/3-dominated regions of NGC 4459 
and NGC 4526 that correspond to the strongest absorption features 
in the unsharp-masked images. 

Regular d ust distrib u tions trace well regular velocity fields, 
consistent with IHo et alJ i20()2l) . All galaxies with perfect dusty 
disks have indeed very regular gas kinematics, with NGC 3032 the 
only possible exception. On the other hand, the absence of regu- 
lar dust lanes does not imply irregular kinematics, as very coher- 
ent gas motions correspond also to less defined dust morphologies 
(e.g., NGC 2974, NGC 4278), and even to quite complex ones (e.g., 
NGC 4150, NGC 5846). 

The unsharp-masked maps also highlight the presence of sep- 
arate stellar components, such as nuclear disks (e.g., NGC 821, 
NGC 5308), peanut-shaped bulges (NGC 2549), and bilobate struc- 
tures reminiscent of a bar (NGC 2699, NGC 4262). Interestingly, 
nuclear disks are a common feature in galaxies with no or only 
weak emission. 



5 KINEMATIC MISALIGNMENT BETWEEN GAS AND 
STARS 

The distribution and kinematics of the ionised gas in early-type 
galaxies have long been known to be often decoupled from that 
of the stars (see Bertola & Corsini 1999, for a review). Figure l4bl 
includes the stellar velocity fields from Paper III to facilitate the 
comparison between the [O III] and stellar kinematics, and shows 
that our sample galaxies are not an exception in this respect. The 
velocity maps show that the motion of gas is often decoupled from 
that of the stars, and that the angular momenta of the gas and stars 
have the same orientation in only a few cases. 

A number of studies have used the dist ribution of the misa lign- 
ment between the spatia l distributions (e.g.. | M artel et all2004l and 
angular momenta (e.g..l KannaDPan & Fabrican ?^20^^ of gas and 
stars to constrain the origin of the ionised gas in early-type galax- 
ies. The orientation of the dust relative to that of the stars has also 
serve d this purpose (e.g.. Ivan Dokkum & Fran?3ll995t iTran et alj 
l200lh . With 48 representative E and SO galaxies in both cluster and 
field environments, the SAURON sample constitutes a good basis to 
test simple hypotheses on the origin of the gas. 

We used the velocity maps of Figure l4atl4bl to accurately mea- 
sure the direction of the stellar and gaseous rotation, and compared 
the distribution for the kinematic misalignments between gas and 
stars with what is expected if the gas has either an internal or exter- 
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Notes: (1) NGC number. (2) Median kinematic misalignment in degrees. 
(3) Standard deviation for the kinematic misalignment in degrees. (4) Max- 
imum distance from the center used to derived the values in Columns (2) 
and (3). For NGC 4150 we excluded also the central 2". 

nal origin. We trace the direction of maximum rotation as a func- 
tion of radius, using the harmonic-ex pansion method for analysing 
two-dimensional kinematics maps of Krajnovic ( 2004). We divided 
the velocity field in concentric circular annuli and fitted the first 
four terms of a Fourier series to the velocity angular profiles in 
each of the annuli. The angular phase of the first order term cor- 
responds to the direction of maximum rotation, ^gas or <^star. The 
width of the annuli increases geometrically to follow the increase 
in the SAURON bin size. The fit automatica lly stops toward t he edge 
of the field or of the gas distribution (see lKrainovidE004l) . In the 
lower panels of Figure l4bl we plot the difference between (^gas and 
<^star as a function of radius. The gas-star kinematic misalignment 
is always shown between ±180°, and is derived only in objects 
with clear gas detection and where stellar rotation is observed. We 
further excluded galaxies without sufficiently extended emission 
(NGC 2699 and NGC 5831). In NGC 4570 and the external regions 
of NGC 7457 the emission is too fragmented for the harmonic fit 
to converge. The complex stellar kinematics in NGC 4550 are also 
hard to describe with this method. 

From the misalignment profiles of Figure l4bl we derived for 
each object a median value and the standard deviation, which are 
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Figure 6. Distribution of tlie median values for tlie kinematic misalignments 
between star and gas for all galaxies listed in Table|2](top) or only for lentic- 
ular (middle) and elliptical (bottom) galaxies. The values of Tablelllare now 
shown between 0° and 180°. 



tabulated in Tabled We considered NGC 4550 and NGC 4570 as 
perfectly co-rotating systems (see AppendijjB}, and used only the 
central gaseous systems of NGC 474 and NGC 2549, as they are 
more likely to be settled. For the same reason we excluded the fil- 
amentary structures extending from the central disks in NGC 5198 
and NGC 5982. In NGC 3414 we considered only the central re- 
gions because of a strong decoupling in the stellar velocity field. 
Table 121 lists the adopted radial ranges. The values for the stan- 
dard deviation of the kinematic misalignment always account for 
the observed twists in the gas and stellar velocity maps since the 
harmonic-expansion generally measures </>gas and </>star to within a 
few degrees. 

Figure|6|presents the distribution of the average values for the 
kinematic misalignment between gas and stars in our sample, now 
shown between 0° and 1 80° . Its principal feature is a pronounced 
excess of objects with gas in prograde orbits with respect to objects 
with gas in retrograde orbits. Figure|6|shows also that elliptical and 
lenticular galaxies have very similar distributions. Despite the mod- 
est number of objects, a KS-test shows that the two distributions 
are identical at a Icr-level (p = 77%). The distributions of field 
and cluster galaxies are also not significantly different (p = 43%). 
The distribution of kinematic misalignments is also independent of 
the galaxy luminosity. 

To interpret Figure |6| let us assume that the intrinsic shape 
of our sample galaxies is mildly triaxial (Franxetal. 1991). In 
this situation, stable closed orbits for the gas are allowed only in 
two planes: the plane containing the short and intermediate axis 
and the plane containing the long and intermediate axis. When 
gas is acquired from random directions, it will settle more often 
in the long-intermediate plane and the chances of settling in the 
short-intermediate plane will scale with the degree of triaxiality 
iSteiman-Cameron & Durisenll982h . 

If the origin of the gas is external, and assuming that pro- 
grade and retrograde gas settle in the same way, the distribution of 



Figure 7. Same as FigurelSlbut now showing flatter and rounder galaxies in 
the middle and bottom panel, respectively. 

the gas-star kinematic misalignments will therefore display three 
peaks: two of equal intensity corresponding to counter- and co- 
rotating gas and stars, and a weaker peak for gas in orthogonal ro- 
tation, assuming the stars rotate along the short axis. In a triaxial 
galaxy, however, the stellar rotation axis can lie anywhere in the 
plane containing the short and long axes. Since gas in equilibrium 
can rotate only around these two axes, intermediate values for the 
kinematic misalignment will also be observed. Projection effects 
will further dilute the distribution of the observed kinematic mis- 
alignments, but overall the resulting distribution will be symmetric 
around 90°, with an equal number of counter- and co-rotating gas 
and stellar systems. 

If the origin of the gas is internal (e.g., from stellar-mass loss) 
the gas will rotate in the same sense as the stars, and therefore the 
distribution for the kinematic misalignments will be asymmetric, 
with values mostly between 0° and 90° . 

The observed distribution in Figure|5|is inconsistent with the 
prediction of either of these simple scenarios. Half of the objects 
show a kinematic decoupling that implies an external origin for the 
gas, but the number of objects consistent with co-rotating gas and 
stars exceeds by far the number of cases with counter-rotating sys- 
tems, suggesting that internal production of gas has to be important. 

We note that our objects form a representative, but incomplete, 
sample of the local early-type galaxies population. The degree of 
incompleteness is known as our targets were drawn from a com- 
plete sample (Paper II). Since galaxies with prograde and retro- 
grade gaseous systems form subsamples that are similarly repre- 
sentative of the local galaxy population, incompleteness does not 
have a significant impact on the previous discussion. Incomplete- 
ness is not responsible for the observed asymmetry in Figure|6| 



5.1 The dependency on galaxy flattening and rotational 
support 

The distribution of the average values for the kinematic misalign- 
ment between stars and gas does not depend on Hubble type, galac- 
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tic environment, or galaxy luminosity. It does, however, strongly 
depend on the apparent large-scale flattening of galaxies. Figure^ 
shows that the roundest objects in our sample (625 ^ 0.2) present a 
more symmetric distribution of kinematic misalignments than flat- 
ter galaxies, which instead host predominantly co-rotating stellar 
and gaseous systems. The two distributions are different at a 99% 
confidence level and there is a 53% probability that the distribu- 
tion of rounder object was drawn from a uniform distribution. The 
galaxies in the two subsamples have no significantly different lu- 
minosities or Hubble types and do not live in different environ- 
ments. The two samples are also equally incomplete. The flat ob- 
jects showing co-rotating gas and stars are also not significantly 
different than the rest of the galaxies in this subsample. 

Since for random orientations fairly round galaxies are likely 
to be almost spherical and hence supported by dynamical pressure, 
rather than by rotation, the degree of rotational support could also 
be important to explain the observed dependency on the galaxy flat- 
tening. A first distinction between slowly-rotating, rounder galax- 
ies and fast-rotating flatter objects, can be made using the classi- 
cal V/o — e diagram recently revised by Binney (2005). This dia- 
gram alone, however, would fail to separate galaxies that are char- 
acterised by overall rotation from non-rotating objects with a fast 
rotating core, because the derivation of V/a includes only a lu- 
minosity weighting. In Emsellem et al. (in preparation) we assess 
the level of rotation support in a more quantitative way, adopting 
a quantity that is closely related to the specific angular momentum 
of a galaxy, thus overcoming the limitations of the V/a — e dia- 
gram. Figure |8| shows the distribution of kinematic misalignments 
between gas and stars in fast and slowly rotating galaxies according 
to the criterion of Emsellem et al. (see also McDermid et al. 2005, 
for an illustration of this two kind of objects). Consistent with our 
expectations, the two distributions are remarkably different (at a 2a 
level, p = 95%), as in the case of flat and round objects (Figure.Q. 
In the context of our simple first-order assumptions these results 
suggest that external accretion of gaseous material is less impor- 
tant than internal production of gas in flat and fast rotating galax- 
ies. On the other hand, the more uniform distribution of kinematic 
misalignments in rounder and slowly rotating objects suggests that 
these objects aquire their gas more often. However, it must be bear 
in mind that the interpretation of the misalignment distribution of 
rounder galaxies can be complicated by projection effects, because 
such objects often host kinematically decoupled core, and by the 
fact that gas is subject to only weak gravitational torques in almost 
spherical objects. 



6 DISCUSSION AND CONCLUSIONS 

We have measured the ionised-gas fluxes and kinematics in 48 el- 
liptical and lenticular galaxies both in cluster and field environ- 
ments, using a novel technique to measure emission lines in galac- 
tic spectra where the stellar and ionised-gas contributions to the 
spectrum are simultaneously described. Extensive simulations were 
performed to test this procedure and assess the detection limits in 
measuring gas emission. 

The excellent quality of the SAURON data and the ability of 
our new method to exploit the entire spectral range allowed us to 
detect emission lines down to an equivalent width of 0. 1 A, which is 
set by the current limitations in describing the spectra of early-type 
galaxies with synthetic and real stellar templates. Due to these limi- 
tations neither the H/3 nor the [N II] lines could be always measured 
without imposing on them the kinematics of the [O III] doublet. In 









■ 


— h- 


-H 1 1 1 1 1 I I 1 1 — 1 1 1 1 1 1 1 1 

fast rotating - 

. 1 . . 1 .1. , '- 




slowly rotating - 



30 60 90 120 150 180 

Figure 8. Same as Figurel^but now showing fast and slowly rotating galax- 
ies in the middle and bottom panel, respectively. 



the case of H/3, independent fits lead to biased gas kinematics that 
are easily recognised across the field of a galaxy. This allowed us to 
identify a few galaxies where the H/3 and [O III] kinematics could 
be independently derived and compared. 

Across our sample, Vgas and (jgas are on average accurate to 
within 14 kms^^and 20 kms~^, respectively. Errors on the flux of 
the [O III] and H/3 lines are on average 10% and 20%, respectively. 
Although the H/3 and [O III] kinematics can be different, imposing 
the [O III] kinematics on the H/3 lines does not dramatically affect 
our ability to measure the H/3 fluxes. This is observed both in the 
simulations and in the objects where independently derived H/3 and 
[O III] kinematics could be compared. On the other hand, relying 
on the detection of [O III] emission before measuring H/3 does limit 
our ability to detect weak emission from H Il-regions, where the 
[O III] lines are dimmer than H/3. 

Emission is clearly detected in 36/48 of our sample galaxies 
(75%) and only 5 objects do not show any significant emission. The 
remaining 7 galaxies exhibit weak [O III] lines only or fragmented 
traces of emission. Among clear detections, a mild dependence on 
the Hubble type and galactic environment is observed, with higher 
detection rates in lenticular galaxies and field objects. Emission is 
found in 20/24 lenticular galaxies in our sample (83%) and in 16/24 
(66%) of the ellipticals. This i s remarkably clo se to the detections 
rates of the imaging survey of lMacchetto et all fl996^. who found 
Ha-l-[N II] emission in 85% of SO and 68% of E. The dependence 
on the galactic environment is similarly marginal, although when 
only the Virgo cluster is considered the fraction of galaxies with 
clearly detected emission drops to 55% (10/18), with just 3/9 ellip- 
ticals exhibiting emission lines. The se 3 objects are also the bright- 
est that we observed in this cluster, ^auer et all fc005') also find a 
significantly lower incidence (33%) of galactic dust in Virgo than 
in the rest of the local elliptical galaxy population (47%). They also 
found dust only in the brightest objects of this cluster. 

The observed emission comes with a variety of resolved distri- 
butions, kinematic behaviours, and [O IIl]/H/3 line ratios. It is very 
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often, although not always, associated with dust. Two interesting 
classes of objects can be recognised. 

The first group show settled gaseous systems where star- 
formation is almost certainly occurring, particularly in circumnu- 
clear regions. The defining properties of this class, which includes 
the SO galaxies NGC 524, NGC 3032, NGC 4459, NGC 4526 
and NGC 5838, are a regular disk-like gas distribution and kine- 
matics, very regular and circularly symmetric dust lanes, and the 
lowest [O IIl]/H/3 ratios in our sample. The kinematic signature 
of emission from H Il-regions in circular motion on the equato- 
rial plane is observed in the independently derived kinematics of 
the H/3 lines, which show faster rotation and smaller velocity dis- 
persions than the [O III] kinematics. The detection in three galax- 
ies of CO emission from dens e molecular clouds (in NGC 3032, 
NGC 4459 and NGC 4526 bv lSaee & Wrobeill989h further sug- 
gests star-formation activity. On the other hand, it is likely that gas 
clouds departing from simple rotation end up getting shocked, in- 
creasing the gas ionisation and the [O III] emission. 

The second group of galaxies, which includes NGC 474, 
NGC 3377, NGC 4262, NGC 4278, and NGC 4546, is charac- 
terised by an integral-sign pattern in the ionised-gas distribution 
and by noticeable twists in the gas velocity maps. In addition, 
NGC 4262, NGC 4278, and NGC 4546 display peculiar asymme- 
tries in their cTgas and [O IIl]/H/3 maps. Since all these objects show 
misaligned stellar photometric and kinematic axes (although only 
mildly for NGC 3377, see Paper I), the observed twisting in the 
ionised gas distribution and kinematics is more likely tracing re- 
gions where gas accumulates while flowing in a non-axisymmetric 
potential rather than a warped configuration. 

The presence of a triaxial structure also appears to be under- 
scored by finding in the outer parts of a galaxy gas emission that is 
misaligned with respect to the main stellar body and that is kine- 
matically decoupled from the gas kinematics in the central region, 
as in NGC 1023, NGC 2549, NGC 7332 and NGC 7457. NGC 1023 
is indeed a well-known barred galaxy (e.g., Debattista et al. 2002), 
NGC 2549 shows a peanut-shape structure in the HST images, 
and the presen ce of a bar in NG C 7332 was extensively discussed 
bvpFalcon-Bar roso et alj j2004h . Taken together with NGC 2974, 
where the gas distribution observed with both SAURON and HST is 
consistent with the presence of nested bars (Krajnovic et al. 2005; 
Emsellem et al. 2003), these patterns show how dramatically the 
gas can respond to the presence of even a modest non-axisymmetric 
perturbation of the gravitational potential. 

Weak bars or unsettled configurations can contribute to ex- 
plain why the ionised-gas kinematics is rarely consistent with sim- 
ple coplanar circular motions. Yet, despite complex structures are 
often observed in the velocity maps, the gas almost never displays 
irregular kinematics and instead generally shows coherent motions 
with smooth variations in angular momentum. In the majority of the 
cases the gas kinematics is decoupled from the stellar kinematics. 

We have measured the kinematic misalignment between stars 
and gas and derived a distribution of mean kinematic misalignment 
values to draw on the origin of the gas. Although half of the objects 
show a kinematic decoupling that implies an external origin for 
the gas, the distribution of the misalignment values between stellar 
and gaseous angular momenta is inconsistent with a purely exter- 
nal origin. In particular, the number of objects with co-rotating gas 
and stars exceeds by far the number of cases with counter-rotating 
systems, suggesting that internal production of gas has to be impor- 
tant. The distribution of the kinematic misalignment between stars 
and gas does not depend on the Hubble type, galactic environment, 
or luminosity of our sample galaxies. It does, however, strongly 



depend on the apparent flattening of galaxies, and their level of ro- 
tational support. 

These results demonstrate that the origin of the gas in early- 
type galaxies is not yet a "solved problem" - more clues are needed. 
Measuring the metallicity of the interstellar medium can provide 
the needed insight. Indeed if the gas originated in stellar-mass loss 
its metallicity should be related to that of the surrounding stars, in 
contrast to what is expected if the gas has an external origin. Yet, to 
measure the chemical composition of the gas it is first necessary to 
understand the ionisation of the gas, i.e. what causes some atomic 
species to emit more efficiently than others. 

In this respect, the SAURON [O IIl]/H/9 maps reveal a wide 
range a values across the sample and considerable structures within 
single galaxies. Despite the limitation of the [O IIl]/H/3 ratio as an 
emission line-diagnostic, this finding suggests either that a variety 
of mechanisms is responsible for the gas excitation in E and SO 
galaxies or that the metallicity of the interstellar material is quite 
heterogeneous from galaxy to galaxy. The [O IIl]/H/3 maps always 
show central gradients, where the line ratios tend to values always 
between 1-3. Since at these [O IIl]/H/9 regimes LINERs and Transi- 
tion objects are the most common class of emission-line nuclei, our 
measurements may be consistent with the finding that the majority 
of E an d SO emission-lin e nuclei belong to these two spectroscopic 
classes JHo et alj[l9973) . In the outer parts of our sample galax- 
ies, except for the objects where star-formation is occurring, the 
[O IIl]/H/3 is also typically between 1-3, although with large scat- 
ter. This could also be consistent with a LINER-like classification 
for the extended emission of early-type galaxies, which has been 
observed in many cases (e.g.. 'Philli ps et al]|l98d : Izeilinger et alj 
1996; Goudfrooii 1999). More emission-line diagnostics such as 
the [N Il]/Ha ratio are needed to confirm the spectral classification 
of the observed emission. 

Excluding a central AGN as the ionising mechanism for the 
emission observed at kilo-parsec scales, shocks, thermal conduc- 
tion, and both young and post-AGB stars are all potential sources 
of ionisation. S hocks are known t o occur while gas flows in barred 
potential (e.g.. lRoberts et ai]|l97 9). therefore their role has to be 
important in galaxies where the presence of non-axisymmetric 
perturbation of the potential is underscored by characteristic pat- 
terns in the gas distribution and kinematics as those discussed be- 
fore. Conductive heating of the warm gas by hot electrons from 
the X-ray emitting gas is a readily available source of energy in 
most early-type gal axies, sufficient to p ower the observed nebu- 
lar emission (e.g.. lMacchetto et alJll99db . High spatial resolution 
X-ray data have also shown a striking coincidence between the 
spatial distribution of the X-ray and ionised-ga s emission (e.g., 
iTrinchieri & Goudfrooiil200llSparks et alj2004l) . supporting a ca- 
sual link between them. Post-AGB stars could represent the most 
common source of ionisation. They can p rovide enough ionising 
photons to power the observe e mission iMacchetto et alJll996h 
and lead to LINER-like line ratios iBinette et alll994lGoudfrooii 
I1999t) . The most compelling piece of evidence supporting this sce- 
nario is the finding that the emission-line flux correlates very well 
with with h ost galaxy s tellar luminosity within the emission-line re- 
gion (.Macchetto et alj|1996i) . This correlation suggests indeed that 
the source of ionising photons are distributed in the same way as 
the stellar population. The SAURON observations further support 
this view, since as noted in i|3.2l the flux distribution of the lines 
closely follows that of the stellar continuum. This is shown by the 
smooth appearance of many EW-maps, in particular for the H/3 
line. Finally, the role of H Il-regions cannot be ruled out only on 
the basis of [O IIl]/H/9 ratios > 1, and should be explored in more 
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details in ligtit of ttie recent claims based on GALEX data of Yi et 
al (2005), that a substantial fraction of nearby early-type galaxies 
recently underwent star-formation activity. 

Complementing this survey with integral-field spectroscopic 
data in the Ha+[N llj wavelength region will mark an important 
step to further understand what powers the nebular emission in 
early-type galaxies, with bearings on the origin of the gas in these 
systems. 
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APPENDIX A: TESTS ON THE EMISSION-LINE 
MEASUREMENTS 

In this appendix, we attempt through a number of experiments to 
establish how accurately the position, width, amplitude, and flux 
of an emission line can be measured as a function of its relative 
strength with respect to both the statistical and systematic devia- 
tions from the fit to the stellar continuum. 



A I The synthetic data 

To simulate the typical situation in the SAURON spectra, we con- 
structed a large number of synthetic spectra, adding to an old stellar 
population varying amounts of H/3 and [O III] emission and statisti- 
cal noise. Specifically we used a 10-Gyr-old solar-metallicity tem- 
plate from the Vazdekis library, to which we added [O III] emission 
lines with A/N ranging from 0.5 to 20, and H/3 emission rang- 
ing from 1/10 to 10 times the [O III] emission. No [N I] emission 
was added. We included Poisson noise to simulate spectra with 
a S/N in the stellar continuum of 200, 100, and 60, the latter 
value being typical of the spectra in the outer parts of our sample 
galaxies. The stellar spectra were all broadened by the same line- 
of-sight velocity distribution, in this case a simple Gaussian with 
a — 150 kms~^centred at zero. The velocity of the [Olll] lines 
ranged from -250 to 250 km s~^, with a constant intrinsic width of 
''"gas ~ 50 kms~^. Since the stellar spectra are kept at their rest 
frame, the gas velocities are effectively relative to the stellar ones. 
We imposed the same kinematics on all emission lines. 

Overall, we created 15390 synthetic spectra, which were anal- 
ysed following exactly the same procedure as for the real data 
( i|2. It . As the stellar kinematics can be measured very accurately at 
the considered S/N and a values, we always adopted the correct v 
and a values while matching the stellar continuum. 



A2 General Results 

Fi gure s IA 1 al A 1 c I show how the measured position, width, ampli- 
tudes, and fluxes of the lines deviate from their input values as a 
function of the measured A/N and for different S/N values. 

The recovery accuracy of the position and width of the lines 
is a strong function of the measured A/N while the amplitudes do 
not seem to depend on it. The quality of the data, as quantified by 
the S/N values, does not seem to affect the Vgas and agas mea- 
surements, while it is dominant in determining how well the am- 
plitudes are estimated. The accuracy in measuring the amplitudes 
of the lines is indeed limited only by the level of statistical noise 
associated to the spectra. 

The precision in estimating the emission-line fluxes turned out 
to be dominated by how well the amplitudes are measured, as it also 
scales mostly with S/N instead of A/N. A dependency on A/N 
is introduced through the impact of the line-widths in the flux mea- 
surements, although only to a lower extent. For the [O III] lines 
this is visible as a weak modulation in the flux recovery accuracy 
with A/N, which indeed does not exactly parallel the small vari- 
ations in the amplitude estimation. The effect on the H/3 line is 
more obvious and is illustrated by a loss of accuracy in measuring 
the fluxes for stronger lines. This is due to synthetic spectra where 
the H/3 line is much stronger than quite weak [O III] lines, so that 
the measured width of the lines is subject to larger fluctuations. 
As the [O III] kinematics is imposed on the H/3 lines while fitting 
them, this can result in an error on the flux measurements that gets 




Figure Ala. Recovery accuracy of the emission-line parameters as a function of the measured A/N for simulations with a statistical S/N = 60. Left panels: 
Difference between the output and input gas velocity (top) and intrinsic velocity dispersion {bottom). The solid and dashed lines indicate the median values 
for the output-input difference and the 68% confidence region around them, respectively. The grey lines show the median values of the formal uncertainties on 
the measured parameters. Middle panels: Same as left, but now for the amplitudes of the [O III] {top) and H/3 {bottom) lines. Right panels: Same as left, but 
now for the fluxes of the [O III] {top) and H/3 {bottom) lines. The accuracy in estimating Vgas and o"gas is a strong function of A/N, while the measurement 
of the amplitudes is not affected by it. The same applies for the fluxes of the [0 III] lines, but for strong H/3 lines these can be subject to larger errors if the 
width of the lines have been estimated from quite weak [O III] lines. In these simulations the main limit to the recovery accuracy of the amplitudes is set by 
the level of statistical noise in the spectra. 



worse the stronger the H/3 emission, which is exactly what is ob- 
served in the simulations. Fortunately, this situation rarely occurs 
in our sample galaxy, since when the [O III] lines are weaker than 
H/3 generally both lines are quite strong. We note that for weak 
H/3 lines the accuracy with which the fluxes are recovered becomes 
again dominated by the error in estimating the line amplitudes. 

Fi gures lAlalAlcl also show how the formal uncertainties on 
the measured parameters closely parallel the standard deviations of 
the measured values from the input ones, demonstrating that our 
error estimates are reliable. Only in the case of the line fluxes, 
our formal error estimates systematically over and underestimate 
by ~ 10% the observed fluctuations for the [O III] and H/3 lines, 
respectively. 

A3 Loss of reliability at low A/N values 

Figures IA 1 al A 1 cl show that as the emission-lines get weaker our 
measurements become less reliable and also start to be dominated 
by systematic effects. 

In the case of the line velocities, there is little evidence for 
a systematic bias, whereas below A/N ~ 2 — 3 the line widths 



start to be overestimated. Two opposite biases are at work in this 
case. As the lines become weaker, their width can either be over- 
estimated as their wings are increasingly lost in the noise level, or 
underestimated as it becomes more likely that the fit just converges 
on pixel-scale random fluctuations. In the latter case the returned 
intrinsic width will always be near zero, which explains the hori- 
zontal stripe of points in the bottom left panels of Fi gures i A I al A iTl 
This last behaviour is seldom observed in real data, as some corre- 
lation between pixels is introduced during the data reduction. As a 
result, the width of the lines tends to be overestimated at low A/N 
regimes. In ilASl it will be shown that template-mismatch can also 
work in this direction. 

Finally, both amplitudes and fluxes tend to be overestimated 
at low A/N regimes. However, except for the fluxes of the [O III] 
lines that can be already biased for A/N ^ 5, this problem is 
not always visible in Figures lAlalAlcl since the measured A/N 
values are themselves overestimated at this level. This is shown 
in the central panels of Fi gure s lAlalAlcl bv the skewed stripes of 
points at very low A/N, which correspond to synthetic spectra with 
lines of the same amplitude. To highlight this bias, in Figure lA2l we 
show the recovery accuracy of the line parameters against the input 
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Figure Alb. Same as Figure lATal but for S/N = 100. Notice how the amphtude and flux of the hnes are less accurately estimated. 



A/N values, clearly not accessible quantities in the real data. In this 
figure both the input-output deviations and the formal uncertainties 
on the amplitudes and fluxes of the lines are shown relative to the 
input values. This illustrates how accurate these measurements can 
be in principle, as even for an input A /N — 3 and a S/N — 60, the 
flux and amplitudes of both lines are recovered to a 30% level, and 
are biased by less than 10%. The H/3 fluxes are less biased than the 
[Olll] fluxes at low A/N regimes (Figure |A2j, as effectively the 
fit to the H/3 line involves a lower number of degrees of freedom. 

When considering the flux and amplitude deviations in lAlal 
lAlcl in relative terms, the observed biases remain quite limited, 
within 10% foi A/N ^ 2.5. 



A4 Assessing the sensitivity to template mismatch: 
simulations 

In real spectra our ability to match the stellar continuum is more 
limited, as the range in metallicities and abundances in our tem- 
plate library do not match what is observed in early-type galaxies. 
This problem can considerably affect our emission-line measure- 
ment, and it is difficult to simulate. Yet, we can understand which 
line is more sensitive to template mismatch by running a second 
set of simulations. Here only the input template is used to match 
the stellar continuum, as opposed to our standard procedure, where 
the whole template library is used. The comparison of these simu- 
lations with the previous set, will tell us if systematics are already 
introduced by allowing the full freedom in our library even in fit- 



ting spectra that can be matched by our templates and to which 
noise and emission lines have been added. 

Figure lA3l shows the result of this experiment. As expected 
the fit to the [O III] lines is not affected by using either the whole 
library of templates or only the true one, while the estimation of 
the H/3 line amplitudes and fluxes is more sensitive in this respect. 
The accuracy in estimating the fluxes of the H/3 lines is on average 
30% worse when allowing the full freedom in the template library, 
with very little dependence on either A/N or S/N, even though 
the models could perfectly account for the stellar continuum. 

When fitting real spectra we can expect that intrinsic limita- 
tions in our templates will make things even worse, although it is 
likely that the impact on the [O III] ineasurements will remain more 
limited. 



AS Assessing the sensitivity to template mismatch: 
[O III]A4959 and [O III]A5007 independent fits 

The [O III]AA4959, 5007 doublet provides a natural testbed to as- 
sess the impact of template-mismatch on our emission-line mea- 
surements using real data. When fitted independently, the [O III] 
lines must indeed return the same velocity and velocity dispersion, 
and their rela tive fluxes sho uld be consistent with theoretical pre- 
dictions fe.g.. lStorev & Zei ppen 2000). 

The previous simulations have demonstrated that the accuracy 
with which the position and the width of the lines are recovered 
scales as A/N^^ . The weaker [O III]A4959 line will therefore dis- 




Figure Ale. Same as Figure lAlal but for S/N = 200 



play much larger deviations from the real gas velocity and velocity 
dispersion than the [O III]A5007 hne. The [O III]A4959 deviations 
will dominate the observed scatter between the velocities and ve- 
locity dispersions derived from each of the [O III] lines. Hence, 
as long as the [O III]A5007 line is sufficiently strong, we can use 
it to trace the gas LOSVD parameters, and use the results of the 
[O III]A4959 fit to study how these parameters are recovered as a 
function of A/ N in this particular spectral region. 

Figure A4 shows the results of this experiment for two rep- 
resentative galaxies in our sample. Plotted against the A /N of the 
[O III]A4959 line are, from top to bottom, the differences between 
the gas velocities and observed velocity dispersions measured from 
the two lines in the [O III] doublet, and their flux ratio. These mea- 
surements correspond only to fits where the A/N of [O III]A5007 
was greater than 5. The most striking feature in Figure IA4| is the 
difference between the behaviours observed in the two galaxies. 
Whereas for NGC 2685 the [O III] lines show position, widths, 
and fluxes that are consistent with our expectations even for very 
weak [O III]A4959 lines, in NGC 4278 the scatter between the 
independently measured LOSVD parameters strongly depends on 
the A/N of the [OlIl]A4959 line. Furthermore, as emission be- 
comes weaker, an increasing fraction of the fits return systemati- 
cally broader [O III]A4959 widths, particularly for y4/7V < 4. 

The observed behaviour is due to template-mismatch. As 
shown by the green lines in Figure lA4l NGC 4278 has much larger 
stellar velocity dispersions than NGC 2685. On the other hand, the 
[O III]A5007 lines have similar widths in both cases. A larger stel- 



lar velocity dispersion a exacerbates the impact of template mis- 
match on the emission-line measurements in two ways. First, more 
massive galaxies have stellar populations with high er metallicities 
and e n hanced abundances of their Q -elements ( e.g..lWor they et all 
| l992t iGreggiol Il997t iTrager etalj r2000: Ma raston et alJ l2003t 
lThoma^^l]l2ob5h . which are more difficult to match with the 
templates in our library. Second, for a given strength of the 
template-mismatch features affecting the emission-line measure- 
ments, in galaxies with larger a the mismatch features are them- 
selves broader and therefore induce a positive bias on the line- 
widths. 



A6 Concluding Remarks 

The previous simulations have shown that the accuracy with which 
the position and width of an emission line can be measured scales 
only with the relative strength of the line with respect to the de- 
viations from the fit to the stellar continuum. On the other hand 
the accuracy in measuring amplitudes and fluxes is dominated by 
the entity of this same residual noise, which itself depends on the 
quality of the spectra and the limitations of the template library. 

The mild dependence of the flux estimates on the precision 
with which the emission-line kinematics is measured is reassuring 
in light of fact that in real galaxies the [O III] and H/3 lines do not 
necessarily share the same kinematics. Fortunately, the loss of ac- 
curacy in measuring the flux of the strongest H/3 lines due to a 




Figure A2. Same as Figure lAlal but for now showing the recovery accuracy against the input A/N values. Furthermore for the amplitudes and the fluxes of 
the lines (middle and left panels), both the input-output deviations and the formal uncertainties are now relative to the input values. 



poorly estimated kinematics from weak [O III] lines is not an issue 
for our sample. 

At low A/N regimes the emission-line measurements become 
dominated by systematics effects and should not be deemed reli- 
able. Even in the framework of ideal simulations the accuracy in 
recovering the H/3 fluxes is hampered by template-mismatch sys- 
tematics. By independently fitting the lines of the [O III] doublet 
we have further investigated in real spectra the impact of template- 
mismatch in the [O III] spectral region. Template-mismatch is more 
important in the presence of larger stellar velocity dispersions and 
can induce to overestimated line- widths. 

For completeness we report that simulations carried out by 
measuring independently the position and width of the H/? and 
[O III] lines show, unsurprisingly, that the H/3 kinematics is subject 
to larger biases than in the case of the [O III] lines. For instance, to 
achieve the same accuracy in recovering the [O III] kinematics for 
A/N = 4, an A/TV = 5 for H/3 is needed. 



APPENDIX B: DESCRIPTION OF INDIVIDUAL 
GALAXIES 

Here we describe the main structures observed in the SAURON 
maps for the ionised-gas distribution and kinematics of each of the 
galaxies in the E/SO sample presented in this paper. We also com- 
ment on the observed dust distribution, the connection to the stellar 
kinematics, the [O IIl]/H/3 ratios, and refer to previous narrow-band 



imaging and long-slit spectroscopic work. We implicitly refer to 
Fieure l4al4bl when describing these structures. 

NGC 474: This well-known shell galaxy ^Tumbul l et ai]ll999l) 

displays a peculiar ionised gas distribution and kinematics. The 
[O III] distribution follows a twisted pattern, along which the [O III] 
lines are stronger at the northern side of the galaxy between 5" 
and 7" whereas they become stronger at the southern end at larger 
radii. Although the H/3 lines are much fainter than the [O III] lines, 
they seem to follow a similar spatial distribution and indeed escape 
detection on the northern side beyond 7". The gas velocity field 
is characterised by a central and an outer component with almost 
opposite angular momenta. In both cases the gas kinematics is de- 
coupled from the stellar one. The width of the lines appear to peak 
in the transition region between the inner and outer gas velocity 
components, which can be explained as a bias introduced by fitting 
a single Gaussian to what in fact in some spectra looks like a su- 
perposition of two separate components with different velocities. 
The absence of strong absorption features is surprising given the 
considerable amount of gas emission. 

NGC 524: This galaxy displays a disk of gas in very regular ro- 
tation, in the same sense of the stars, consistent with the IFS obser- 
vations of Sil'chenko 1 2000). The emission is close to our detection 
limit, however, which explains the patchy appearance of the maps. 
A similarly fragmented distrib ution is shown by the Hcn-[N II] 
narrow-band images of Macchetto et alJ il996h . The gas distribu- 
tion is also matched by an extremely regular distribution of the dust. 




Figure A3. Same as Figure lATal but for models where only the right, input, stellar template have been used to match the stellar continuum. Notice how the 
accuracy in any of derived parameters for the [O III] lines has not substantially increased with respect to those shown in Figure lAlal while the improvement is 
cleai'ly visible for the H/3 amplitudes and fluxes. 



organised in concentric circular lanes. The strength of the H/3 and 
the [O III] emission are comparable, although at the edge of the 
disk the [O III] lines become stronger. 

NGC 821: No significant emission is detected in this close to 
edge-on galaxy. The HST unsharp-masking i mage reveals a nu- 
clear d isk structure but no dust, consistent with lRavindranath et alj 
<200ll) . 

NGC 1023: This galaxy displays [O III] emission that remains 
strong in the outer regions, with a distribution that appears to be 
skewed with respect the galaxy major axis. Yet, the gas shows a co- 
herent, though not perfectly regular, velocity field with a reasonably 
defined zero-velocity curve, which indicates a mild kinematic mis- 
alignment between gas and stars. The approaching south-eastern 
arm indicates that the angular momentum of the gas may vary in 
the outermost parts of the galaxy. This would be consistent with 
th e complex neutral h ydrogen morphology and kinematics found 
bv'Sanci si et alj|l984) . suggesting an interaction with three nearby 
companion galaxies. No noticeable absorption feature is visible in 
the HST image, which also in this case reveals the presence of a 
nuclear stellar disk. 

NGC 2549: Another object with strong and misaligned emission 
in its outer parts, where the gas kinematics is misaligned with re- 
spect to the stellar kinematics. Towards the centre, however, the an- 
gular momentum of the gas quickly aligns itself behind the stellar 
angular momentum and along the galaxy minor axis, suggesting the 



presence of a gaseous system well settled on the galaxy equatorial 
plane. Also in this case no dust is observed, although the unsharp- 
masked image reveals the clear peanut-shape signature of a boxy 
bulge. 

NGC 2685: Also known as the Helix galaxy 
CBurbidae & Burbidej ll959h . this famous peculiar object shows 
a gas distribution almost perpendicular to the galaxy major axis. 
A spiral-like structure is evident in the flux maps, in particular 
on the north-eastern side of the galaxy. This feature could also 
be regarded as a warped configuration where gas is settling on 
the galaxy equatorial plane towards the central regions. The gas 
kinematics is also strongly misaligned with respect to the stellar 
kinematics, on average by 75°, and shows a complex velocity 
field, with a zero-velocity curve that is skewed with respect to the 
direction along which the most extreme velocities are observed. 
Figure I4bl also shows that towards the centre the gas angular 
momentum is aligning itself behind the stellar one. Near IR 
narrow-band HST images suggest the presence of a nuc lear Pa 
disk aligned with the galaxy major axis jBoker et alll999l) . which 
would confirm the presence of a settled central gaseous system. 
The velocity dispersion map also shows peculiar features ±5" 
away from the centre along the minor axis, where cTgas increases 
up to 150 kms~^. No strong evidence for complex line profiles 
is found. The unsharp-masked image reveal a nuclear disk and 
almost perpendicular dust lanes both across the centre and towards 
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Figure A4. Comparison between emission-line parameters measured independently from the [O in]A4959 and [O III]A5007 lines, as a function of the A/N 
of [O III] A4959 and for the case of NGC 2865 (left) and NGC 4278 (right). The top and middle panels show the difference between the measured gas velocities 
and observed velocity dispersion. The lower panel shows the ratio between the fluxes of the [O III] A4959 and [O III] A5007 lines. In each panel the horizontal 
dashed line indicates our expectations (i.e. no difference between the velocities and velocity dispersions and a flux ratio of 0.33), while the solid and dashed 
red lines show the median values for the plotted differences or ratios and the 68% confidence region around them, respectively. In addition, in the middle 
panels the blue dotted line indicates the gas velocity dispersion as traced by the median width of the [O III]A5007 lines, while the green dashed line shows the 
median stellar velocity dispersion at the location of the emission-line fits. 



the north-eastern side of the galaxy, where also the H/3 emission 
is stronger than the [Olll]. Elsewhere the distribution of the 
[O IIl]/H/3 ratio is complex. 

NGC 2695: No emission is detected in this galaxy, and no HST 
image is available. The SAURON unsharp-masked image does not 
show peculiar features, besides highlighting a twisting of the inner 
isophotes. 

NGC 2699: This object shows emission only in the central re- 
gions, making it difficult to judge the gas kinematics. The HST 
image reveals a bilobate structure underscoring the presence of a 
nuclea r bar, and a very small nuclear dust disk that was noticed 
also bv lTranetal]j200ll) . 

NGC 2768: This gal axy is know n for hosting a central dust lane 
along the minor axis jKimlll989h and for showing a kinematic 
decoupling between stars and gas jBertolaetal]|l99A that was 
later interpreted as the result of a gaseous polar-disk structure 
jpried & Illingworth f'l994y The SAURON maps further reveal re- 
markably different distributions for the [O III] and H/3 emission. 
The [O III] lines follow a twisted pattern that is misaligned by 
~45° from the direction of the dust lanes, while the H/3 distribution 
appear s to be rounder HST narrow-band images by Martel et al. 
J^^^ show that the Ha-l-[N II] emission in the central arcsecond 
also follows an integral-sign pattern that is decoupled from the ori- 
entation of the dust. The [O III] velocity field is well defined, show- 
ing peculiar and almost symmetric twists at ~10" from the cen- 
tre. Overall, the gas motions are almost perpendicular to the stellar 
ones. The cTgas-map shows an elongated peak parall el to the zero- 
veloci ty curve, consistent with the observations of iBertola et alJ 
il995l) . and possibly also a rise of cTgas corresponding to the pe- 
culiar twists in the velocity field. The crisscross appearance of the 



[O IIl]/H/3 maps can be interpreted in light of the different shape of 
[O III] and H/3 distributions. 

NGC 2974: SAURON obs ervations for th i s gala xy were previ- 
ously discussed in detail in lKrainovic et alj i2005l) . Here we note 
how the iSVF-maps highlights the presence of a nuclear bar sur- 
rounded by a ring and of two spiral arms connecting to the 
outer parts of the field where emission remains strong. The pres- 
ence of nested bars within this galaxy is discussed in the light 
of the SAURON and H ST observations by Krajnovic et al. and 
lEmsellem et alj i2003t) . The inner bar is responsible in particu- 
lar for the observed twist within the central 4" in the velocity 
field, which is otherwise very regular. The velocity dispersion 
of the gas rises up to 250 kms~^ and remains moderately high 
(~100 kms~^) also in the rest of the map. The possibility of a 
dynamical supp ort for th e gas motions was first investigated by 
^inzano & van der Mareil il994) and recently also by Krajnovic et 
al. Although not perfectly regular, the dust distribution appears to 
trace the general orientation of the gaseous disk, following in par- 
ticular part of the EW-vmg. 

NGC 3032: Although heavy binning limits our ability to judge 
the regularity of the gas distribution and kinematics in this galaxy, 
overall the gas appears to be consistent with a disk of gas in circular 
rotation, in almost the opposite sense than the stars. The presence 
of very regular dust-lanes and of a more regular H/3 kinematics 
(Figure IC let further suggest a disk distribution for the gas. This 
dusty galaxy features the strongest H/3 emission in our sample, and 
the lowest values for the [O IIl]/H/3 ratios. The relative strength of 
the H/3 line compared to the [O III]A5007 line shows a clear ra- 
dial pattern: from twice as bright at the edge of the disk, to ten 
times stronger at ~4", and finally giving way to the [O III] lines 
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very sharply in the cen tral 2". The detection of molecular gas by 
ISase & Wrobell ^98^ further suggest ongoing star formation ac- 
tivity. 

NGC 3156: This dusty galaxy shows a complex distribution for 
the [O III] emission, with filamentary structures extending perpen- 
dicularly to the galaxy major axis. On the other hand, the H/3 emis- 
sion follows more closely the stellar distribution and forms a nu- 
clear ring corresponding to the location of the central absorption 
features visible in the HST images. Despite the peculiar features in 
the [O III] distribution the gas kinematics show rather coherent mo- 
tions and very small line widths. The [O IIl]/H/3 map displays very 
large fluctuations, with high [O IIl]/H/3 ratios 3) corresponding 
to the main [O III] structures and lower values corresponding to the 
nuclear H/3 ring. 

NGC 3377: SAURON observations for the ionised-gas distribu- 
tion and kinematics of NGC 3377 were previously presented in Pa- 
per I. Here we note how the ionised-gas distribution displays an 
integral-sign shaped pattern that is parallelled by a similar twist in 
the gas velocity field. Along this pattern the emission-line fluxes 
display a much milder decrease with galactocentric distance than 
the stellar continuum, resulting in decreasing EW values towards 
the centre where the weaker H/3 emission eventually escapes de- 
tection. The single dust filame nts visible on the n orthern side of the 
galaxy was also recog nised bv lTran et alJ J200lh . 

NGC 3379: The largest member of the Leo Group, this galaxy 
displays a central disk of gas and a number of isolated sources of 
emission. The direction of maximum rotation of the central disk 
is consistent with the orientation of the nuclear dusty ring visible 
in the HST images, and indicates a ~45° kinematic misalignment 
between gas and stars. The impact of template-mismatch in bias- 
ing the (jgas measurements towards overestimated values ( i|3.3.2t 
is visible at the edge of the central disk, where the strength of the 
emission lines approaches the detection limit. On the other hand, 
all isolated sources show very narrow [O III] lines and undetected 
H/3 emission, consistent with large [O IIl]/H/3 ratios. These charac- 
teristics suggest that these sources are in fact PNe. The two sources 
with the largest EW values ~20" south and south-west from the 
centre, h ave positions and velo cities consistent with the PNe #4 
and #7 of lCiardullo et alj jl993l) . respectively. 

NGC 3384: This galaxy shows weak [O ill] and H/3 emission 
close to our detection limit. The unsharp-mas ked image suggests 
the pr esence of two, nested, disk structures. iRavindranath et all 
yOOl|) noticed the same features in near-IR HST images. 

NGC 3414: This peculiar SO galaxy (Arp 162) shows a spiral 
pattern in the distribution of the gas, which is more evident in 
the EW maps and is accompanied by a 90° twist of the veloc- 
ity field. In the central regions the gas rotates almost perpendicu- 
larly w ith respect to the stars, as noticed by Sil'chenko & Afanasiev 
y004|). The (Jgas -map show a very sharp peak, which is elongated 
in the same direction of the zero-velocity curve, as in the case of 
NGC 2768. Weak filamentary absorption fe atures are visible in the 
southern part of this galaxy, consistent with lTran et aljiJoOll) . 

NGC 3489: Similarly to NGC 3156, also this dusty galaxy fea- 
tures a complex [O III] distribution with filamentary structures and, 
despite that, a fairly regular velocity field and very narrow lines. It 
also shows a more uniform distribution for the H/3 emission, which 
explains the resemblance between the maps for the [O IIl]/H/3 ratio 
and the EW of the [O III] A5007 line. 

NGC 3608: As in NGC 3414, this galaxy shows a spiral pattern 
in the ionised-gas distribution. In this case, however, the emission 
lines are much weaker and no gas is detected between the arms. 
Furthermore both NGC 3414 and NGC 3608 show little overall 



stellar rotation except in the central regions, where the kinematics 
of gas and stars are strongly decoupled. No clear absorption fea- 
tures are visible in the HST image, although iTomita et alJ i200Cl) 
claim to detect a dusty disk. 

NGC 4150: This object displays a extremely smooth EW-maps, 
exemplifying how closely the gas emission can follow the stellar 
continuum in early-type galaxies. Despite heavy binning the Vgas- 
and CTgas-maps suggest fairly regular gas motions in a dynamically 
cold disk. A complex dust structure in the central 3"corresponds 
to the region where the stellar kinematics reveal counter-rotating 
structure (Paper III) and to a minimum in the [O IIl]/H/3 distribu- 
tion. Outside these central features, the gas and stellar kinematics 
are mildly misaligned by 20° . 

NGC 4262: Like NGC 3377, this strongly barred galaxy shows 
an integral-sign pattern in the gas distribution and a twisted V^as- 
field. In addition, this object shows peculiar asymmetric distribu- 
tions for the central values of agas and for the [O IIl]/H/3 ratio. The 
gas and stellar kinematics are strongly decoupled. The unsharp- 
masked SAURON image shows a bilobated structure reminiscent of 
that observed at smaller scales in NGC 2699. 

NGC 4270: This galaxy shows only fragmented traces of emis- 
sion and no evidence for significant absorption features. 

NGC 4278: This galaxy shows the strongest emission, in terms 
of the EW of the lines, in our sample, allowing to appreciate in 
detail the peculiar distribution and kinematics of the gas. As in pre- 
vious objects, the ionised gas displays an integral-sign pattern, a 
twisted Igas -field, and an asymmetric (jgas central peak and distri- 
bution of the [O IIlJ/H/3 ratio. Dust absorption features are visible 
only on the north-western side of the galaxy. Their morphology 
suggests an orientation similar to that of the integral-sign pattern. 
The gaseous and stellar kinematics are misaligned by increasingly 
wider angles, as the stellar and gaseous velocity fields twist in op- 
posite direction towards the outer parts of the fi eld. This object also 
hosts an outer H I ring teaimond et al."'l98l') that may be physi- 
cally connected with the ionised gas give n the sim ilar orientation 
of the two gaseous systems <Goudfrooii et aljl994h . 

NGC 4374: K nown for its low-luminosity radio jet (e.g., 
iLaing &. Bri(U3 IT987^. also this galaxy shows intense ionised-gas 
emission. This is confined along a lane running across the galaxy 
from east to west, which broadens up in the north-south direction at 
large radii. The central dust lane also runs across the nucleus in the 
same direction of the ionised-gas, and perpendicular to the radio 
jet. The gas kinematics indicate a coherent rotation, except in the 
eastern direction where the initially approaching gas gradually re- 
verse its angular momentum and eventually, beyond ~15" is fully 
receding. Clear stellar rotation is observed only in the central ~3", 
in almost opposite direction with respect to the gas. The [O IIl]/H/3- 
map shows rather uniform values, between 1-2. We note that Chan- 
dra observati ons reveal an H-sh aped distribution for the soft X- 
ray emission iFinoguenov & Jones.2001.) that is broadly consistent 
with the ionised-gas distribution. 

NGC 4382: Only weak [O ill] emission is detected in this galaxy, 
mostly towards the centre. The HST unsharp-masked image sug- 
gests the presence of a dust lane r unning across the n ucleus along 
the minor axis direction. However. lLauer et alJi200^ show the ab- 
sence of dust reddening in the centre and find that NGC 4382 dis- 
plays a central minimum in its surface brightness distribution. 

NGC 4387: This galaxy shows only marginally detectable emis- 
sion in few spectra, corresponding to isolated sources, possibly 
PNe. 

NGC 4458: This galaxy displays only fragmented traces of emis- 
sion, in particular close to the centre. The HST unsharp-image re- 
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veals a nuclear disk jMorelli et all2004) . consistent with the pres- 
ence of a central fast-rotating component in the stellar velocity 
field. 

NGC 4459: This SO galaxy shows an extremely well-defined 
dust distribution and regular gas kinematics. The H/3 emission 
peaks in a circumnuclear ring, corresponding to the region where 
the dust distribution appears the most circularly symmetric and 
flocculent. On the other hand, the [O III] emission very closely fol- 
lows the stellar continuum. The values of [O IIl]/H/3 ratio therefore 
show a very similar radial pattern to NGC 3032, with a minimum 
in a circumnuclear region and stronger [O IIl]/H/3 ratios toward the 
centre and the edge of the gas disk. The independently derived H/3 
kinematics show faster rotation and smaller velocity dispersions 
than the [O III] kinematics, suggesting that the H/3 line is tracing 
better than the [O III] doublet a kinematically colder component of 
gas, which is currently forming stars. The detection of molecular 
gas by Sage & Wrobel ( 1989) supports this picture. 

NGC 4473: Also this galaxy displays scattered traces of emis- 
sion, mos tly from the [Ol l l] dou blet. [N II]+ Hq emission was 
found by iMacchetto et alJ h996h but not by lOoudfrooii et alJ 
il994) . The HST unsharp-masked image highlights the presence 
of a nuclear disk, already noticed by van den Bosch etal, 1 1994) . 
No dust is visible, contrary to the finding s of Ferra ri et alJ i 19991) 
but consistent withl van Dokkum & Franxi tl995iV 

NGC 4477: lonised-gas emission in this object is very concen- 
trated toward the centre whereas it decreases more gently than the 
stellar continuum at outer radii. The gas kinematics is rather regu- 
lar, and misaligned with respect to the stellar kinematics by ~30° . 
H/3 and [O III] emission have comparable strengths over most of 
the field except in the central 5" where the [O IIl]/H/3 ratio rises 
to values ~3. Beside the evident central absorption features, weak 
dust lanes are visible ~7" to the south. 

NGC 4486: The central member of the Virgo cluster, M87 dis- 
plays the most complicated gas structures in our sample. The 
ionised-gas distribution features several filamentary structures and 
a very bright central co mponent, consistent with previous narrow- 
band observations (e.g., ISparksetalj|l993h . In the central 2" the 
emission lines have complex profiles that are poorly matched by 
single Gaussians, which explains the large central values for (Tgas. 
No [O III] emission is detected along the radio jet, whereas we 
found H/3 emission when this line is fitted independently. The 
gas velocities show extreme variations across the field, particularly 
along the main filamentary structures to the east of the nucleus. The 
ionised gas is outflowing from the centre to the north, decelerating 
at larger radii. The main features of the Vgas -map are remarkably 
consistent with the long-slit observations of ISoarks et alj il993h . 
Most of the observed emission comes with rather intermediate val- 
ues for the [O IIl]/H/3 ratio, between 1 and 2, except in the south- 
eastern filament and on the closer end of bright patch ~30"south- 
east from the centre. 

NGC 4526: Another prototypical example of an early-type 
galaxy with a perfect dust disk, very regular gas kinematics, a 
circumnuclear region where H/3 emission is particularly strong, 
and an independently derived H/3 kinematics indicating faster ro- 
tation and smaller velocity dispersions than the [O III] kinematics. 
As in NGC 3032 and NGC 4459 molecular gas was detected by 
ISage&Wrobei l'l989^. 

NGC 4546: Another galaxy showing an integral-shape pattern in 
ionised-gas distribution, along which the emission is stronger on 
the western side of the galaxy ~10" from the centre, and on the 
eastern side beyond 15". As in NGC 4262 and NGC 4278 we ob- 
serve a strong twisting of the gas kinematics, and an asymmetric 



distribution for the central (jgas values and for the [O IIl]/H/3 ra- 
tios. The gaseous and stellar moti ons are strongly decoupled, with 
almost opposite angular momenta <GallettJl987h . Absorption fea- 
tures are found only close to the centre. 
NGC 4550: Known for hosting two counter-rotating stellar disks 
iRubin et alJl992tlRix et aljl992 l). the ionised-gas distribution and 
kinematics of this peculiar object, together with the observed dust 
morphology, is consistent with gas moving in a fairly inclined and 
dynamically cold disk. A spiral arm extending to the south is ev- 
ident in the EW maps. IWiklind & Henkell i200ll) note a similarly 
lopsided distribution for the molecular gas they detect. Consistent 
with the measurements of Rubin et al. and Rix et al., the ionised- 
gas rotates against the stars in the disk dominating the stellar kine- 
matics along the major axis beyond 10". On the other hand, the gas 
rotates in the same sense of the stars in the main body of the galaxy. 
Since this implies that the gas does not necessarily have an external 
origin, we decided to adopt a value of zero for the kinematic mis- 
alignment between gas and stars in NGC 4550 for the purpose of 
the analysis of fjs] Independently derived maps for the [O III] distri- 
bution and kinematics were presented by Afanasiev & Sil'chen^ 
(2002). 

NGC 4552: Ionised-gas emission, mostly from the [O III] dou- 
blet, is clearly detected in the central regions of this object al- 
though template-mismatch severely affects the measured values of 
o-gas. In the central 4" the gas kinematics is strongly decoupled 
from the stellar kinematics, confirming that the observed emission 
is not purely an template-mismatch artifact. Ha+[N II] emission 
was detected also by Macchetto et al. ( 1996^). Nuclear [O III] emis- 
sion was found also by Cappell ari et alJ 11199^. Patchy absorption 
features are visible on the north-western side of the centre and ~8" 
to the north-east. Upon closer inspection, the HST image shows the 
presence of a very small ring of dust (r~0'.'25) surrounding the 
nucleus, consistent with the analysis of Carollo et al. 1 1997). 

NGC 4564: No significant emission is detected in this galaxy. 

NGC 4570: This object displays several scattered sources of 
emission. The observed emission is characterised by narrow lines, 
which appear to trace gas clouds coherently rotating in the same 
sense of the stars. We found no tr aces of dust in the HST ima ge, 
where the nuclear disk studied by Ivan den Bosch et all il998h is 
clearly visible. 

NGC 4621: Emission in this galaxy is found the centre and in a 
number of isolated sources. Also in this case the scattered sources 
broadly follow the main stellar stream. As in NGC 4552 template 
mismatch affects the central o-gas measurements, but in this case the 
gas velocities do not show significant rotation and instead appear to 
be offset from the central stellar velocities, which are close to zero. 
The presence of nuclear emission shou ld therefore be confirme d. 
No Ha+[N II] emission was detected bv lGoudfrooii et al]^1994^ 

NGC 4660: This galaxy shows clearly detected [O III] emission, 
possibly originating from PNe, in a few spectra, in particular from 
the southern part of the gala xy. A nuclear disk i s visible in the HST 
image, as recognised also bv lLauer et alj i2005h 

NGC 5198: The ionised-gas distribution in this object is charac- 
terised by a central component and a filamentary structure extend- 
ing from it to the north. Near the centre the gas rotates rapidly and 
almost perpendicularly with respect to the stars in the kinematically 
decoupled core (Paper III). No coherent pattern for the gas motions 
in the filaments is recognised, except that the gas move s relatively 
slowly . No evidence for dust is found, consistent with lTran et alJ 

NGC 5308: No emission is detected in this galaxy, which is al- 



48 Sarzi et al. 



most certainly seen edge-on given tlie extreme inclination of its 
nuclear stellar disk (Krajnovic & Jaffe 2004) 

NGC 5813: The ionised-gas distribution in this object is confined 
within a lane running along the galaxy minor axis. The observed 
gas kinematics is rather complex outside the central ~ 8", but in the 
central regions the gas show coherent rotation in approximately the 
same sense of the stars. The present V^as-map is consistent with the 
measurements shown in Paper II, once the effect of spatial binning 
is accounted for. Dust features take the form of filaments running 
towards the centre from the south and the east. Absorption is also 
strong on the north-western side of the nucleus. The emission-line 
fluxes have a similar distribution near the centre. In particular the 
flux distribution is lopsided along the major-axis direction. Note- 
worthy are also the strong central o-gas peak and the rather uniform 
distribution of intermediate values for the [O III]/H/? ratio. 

NGC 5831: This galaxy shows emission only close to the centre, 
making it hard to comment on the gas kinematics. No trace of dust 
absorption features is found. 

NGC 5838: Another SO galaxy with a very regular dust disk. 
Considering the small scale of the disk and the effects of atmo- 
spheric blurring, the observed gas distribution and kinematics is 
remarkably similar to that of NGC 4459 and NGC 4526. 

NGC 5845: Only weak [O ill] emission is detected towards the 
centre of this galaxy. The HST imag e shows a nuclear stellar disk 
with an inner dust ring, as noticed bv lTranetaljj200lh . 

NGC 5846: Ionised-gas in this galaxy is found within ~IO"and 
along a filament running to the north-east. Although characterised 
by large fluctuations in some places, the Vgas field clearly indicates 
a nearly cylindrical rotation for the gas, in the opposite sense of the 
stars. Remarkably, the zero-velocity curve runs across the north- 
eastern filament, suggesting this is not a separate component from 
the gas closer to the centre. The HST image shows dust along an 
almost parabolic lane swinging around the nucleus. Upon closer in- 
spection, the emission-line fluxes follow very closely the dust dis- 
tribution, in particular from H/9. Relatively intermediate [O IIl]/H/3 
values are observed across most of the field. 

NGC 5982: As in NGC 5198, the gas in this bright elliptical 
is found in a rapidly rotating central component and in a fila- 
ment extending from it, in this case to the south. In NGC 5982, 
however, the motions along the filament are remarkably coher- 
ent. The gas velocities seem to suggest that the gas, initially ap- 
proaching the observer from the southern end of the filament, 
progressively decelerates and plunges back toward the centre fol- 
lowing a curved trajectory. The gas in the centre rotates in the 
same sens e as the stars in the well-known kine matically decou- 
pled core iWagner. Bender & Moellenhofj|l988h . Absorption fea- 
tures are neither found i n the V-band unsharp -masked image nor in 
the V-H colour image of lOuillen etlj]<2000h . 

NGC 7332: SAURON maps for the ionised-gas distribution and 
kine matics for this SBO galaxy were already presented by 
iFalc on-Barros o et alJj2004h . Here we note how, like in NGC 2549 
and NGC 1023, the gas emission remains relatively strong in the 
outer parts of the galaxy, where the gas distribution is misaligned 
with respect to the galaxy major axis and the gas kinematics is de- 
coupled from that observed in the inner regions. We also note the 
large variations in the [O IIl]/H/3 values, and the presence of dust 
in the HST-WFPC image. 

NGC 7457: Another galaxy with strong and misaligned emission 
in its outer parts. Near the centre the gas rotates in the opposite 
sense to the gas in the outer parts and also against the motion of 
the bulk of the stars in this object. No dust is found in the HST 
unsharp-masked image. 



APPENDIX C: MAPS FOR INDEPENDENTLY DERIVED 
[O III] AND H/3 KINEMATICS 

In this appendix we show for 10 galaxies maps for independently 
derived velocity and intrinsic velocity dispersion of the [O III] and 
H/9 lines. These objects are the only ones in our sample for which, 
over most of the regions where emission is observed, the measured 
H/3 kinematics do not suffer from the template-mismatch effects 
described in i|2.1l This was ensured by comparing the H/3 velocity 
maps to the stellar velocity maps and by further noticing that spu- 
rious H/3 emission lines tend to be considerably broader than the 
[O III] lines. 

Figures IC 1 alC I cl show the independently derived [Olll] and 
H/3 velocity and intrinsic velocity dispersion of these 10 galaxies. 
Since among them, 3 are giant ellipticals and 4 show regular and 
concentric dust lanes, Figure lcTal shows first the maps for the re- 
maining 3 galaxies, while figures ICIbl and ICIcI shows maps for 
these two groups. In addition, Figure lCTal shows again the case of 
NGC 3489, to highlight the main characteristics of an unreliable 
H/3 kinematics. 
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Figure Cla. Maps for independently derived [O III] and H/3 kinematics. From top to bottom: i) and ii) mean velocity of the [O III] and H/3 lines iii) and iv) 
intrinsic velocity dispersion of the [O III] and H/3 lines. The maps show only regions where both [O III] and H/3 lines were detected following our standard 
approach. Left: the case of NGC 3489 is shown again to highlight the main characteristics of an um'eliable H/3 kinematics. In addition to the velocity bias 
discussed in i|2. II the spurious H/3 lines are considerably broader (by ~ 50km s~^) than the [O III] lines. Right: three galaxies with reliable independently 
derived H/3 kinematics over most of the regions where emission is observed. In the central regions of all three galaxies the [O III] lines trace faster rotation 
velocities than the H/3 line. In NGC 3414 and NGC 4278 the central Cgas increase of the H/3 lines appeal's also to be much less dramatic than the o-gas increase 
of the [O III] lines. For NGC 3414, template-mismatch is affecting the H/3 kinematics measured in the outer parts of the maps. 
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Figure Clb. Same as Figure lcTal but now showing giant elliptical galaxies. Notice the naiTower width of the H/3 Hnes at all radii. 
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Figure Clc. Same as Figure IcTal but now showing galaxies with circularly symmetric dust lanes. In all cases the H/3 kinematics show faster rotation and 
smaller velocity dispersions than the [O III] kinematics. 



